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1.0 INTRODUCTION AND SCOPE OF DOCUMENT 


1.1 ASCENT/DESCENT OUTPUT OVERVIEW 

The Ascent/Descent output product is produced in the OPIP routine from three 
files which constitute its input. One of these, OPIP. IN, contains mission 
specific parameters. Meteorological data, such as atmospheric wind 
velocities, temperatures, and density, are obtained from the second file, the 
Corrected METeorological DATA file (METDATA). The third file is the 
TRJATTDATA file which contains the time-tagged state vectors that combine 
trajectory information from the Best Estimate of Trajectory (BET) filter, 
LRBET5 (Applicable Document 4, Section 3.0), and Best Estimate of Attitude 
(REA) derived from IMU telemetry. 

Ascent/Descent output products are provided with two files. The major product 
contains the BET, which is in the file BETDATA. The other product file 
contains the system parameters which were used to determine the BET. This is 
The Navigation Block, or NAVBLK, file. 

Each of the product files is delivered on three magnetic computer tapes. Two 
of these are in binary format, one of which is formatted to be read by a CDC 
Cyber machine, and the other is formatted to be read by a Sperry Uni vac 
machine. The third tape is written in ASCII format and is used to produce a 
microfiche display of the output. 


1.2 SCOPE OF DOCUMENT 


This engineering description defines each term in the two output data files. 
The description of the BETDATA file includes an outline of the algorithm used 
to calculate each term. Most computations are performed in the program OPIP 
or in simple subroutines called by OPIP. Some subroutines are extensive, 
however, and their algorithms are described separately in a section following 
that for OPIP. 
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To facilitate describing the algorithms, a nomenclature is defined in Sections 
2.1 through 2.4. The description of the nomenclature includes a definition of 
the coordinate systems used. 

The NAVBLK file contains navigation input parameters. Each term in NAVBLK is 
defined, and its source (e.g., user input) is listed. The production of 
NAVBLK requires only two computational algorithms. These two algorithms, 
which compute the terms DELTA and RSUBO, are described in Section 4.2.2 
following the definition of NAVBLK terms. Finally, the distribution of data 
in the NAVBLK records is listed. 

The English system of units is used throughout this document. Unless 
specified otherwise, lengths, velocities, and accelerations are output in 
feet, feet/sec, and feet/sec^, respectively. Mass is output in slugs, force 
(and weight) in pounds, temperature in degrees Rankine, and angular measure 
(eg., latitude, longitude, Euler angles, angle of attack and sideslip angle) 
are in degrees. Also, unless specified otherwise, time is output in seconds. 
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2.0 NOMENCLATURE 


2.1 VECTORS 




will represent a three-element position vector 
from point A to point B; i.e., the position of 
B with respect to A. The components of the 
vector are with respect to the reference frame 

A. A. 

F whose orthogonal basis set is [ip, jp , kp]; 

A A. 

ip, jp, and kp are unit vectors. Thus, 


ki 


rA[F] 

B 


Ta[F] - Xgip + Ygjp + Zgkp. 


?a[f] 


• f'/\CF] is the unit vector in the direction of r^[F]; r^[F] 

B B B 


I?aCf]| 

B 


• ^ ^ which is a velocity if F is an inertial frame. 


B 


B 


2 

• r.[F] = ^— 2 r r.[F] which is an acceleration if F is an inertial frame. 

B B 


• Vp[F] If a relative velocity (say, wind relative to some point, P) is 
wind 

indicated or if it is not useful to express a velocity in terms of 

a time derivative, then the velocity may be expressed as Vp[F] 

wind 

where, again, F is the coordinate frame in which the vector 
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components are expressed. 


• Likewise, the acceleration of a point P with respect to a point A 

p 

as expressed in frame F is as shown. In this document, the term 
"contact acceleration" (which is the term most often used in the 
software listings) is equivalent to sensed acceleration. 


• §[M50] is the M5D state vector. 


• o(w) will represent the standard deviation of the components of the 

vector w. The term "uncertainty" will be equivalent to the standard 
deviation. If w has three components, then so will S(w). That is, 
if w = w^i + w^j + w^k, then a(w) = [a(w^), o(Wy), o(w^)]^. 


• w[F] is the angular velocity expressed in the F-frame. 


• i^[F] is the angular acceleration expressed in the F-frame. 


2.2 MATRICES 

# [FI ->■ F2] represents a transformation matrix which transforms a vector 
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from frame FI to frame F2. Thus, 


= [FI - F2] f^[Fl] 

B B 

In terms of elements, this matrix may be expressed as 
[FI F2] = M = [m. .] where i represents the row and j 

’ \j 

represents the column of element m^j. 


• [FI > F2f represents the transpose of the above matrix. Thus, 

since 

[FI > F2]^ = [F2 > FI] , 
then 

= [FI » F2f f^[F2]. 

B B 

• (|:(w) represents the covariance matrix of the vector w. If w has 

n elements, then (f(w) will be an n x n matrix. 


2.3 SPECIFIC LOCATIONS 

The following convention will be used to express the following points, or 
locations. 

• CM = the center of mass (usually equivalent to the center of gravity). 
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• NB = the Navigation Base (= Nav Base). 

• (9 = the center of the Earth. 

• L = origin of Landing Field (Runway) coordinate system. 

• S = origin of either Launch Site for Ascent analysis or Landing Site 

Runway threshold for Descent analysis. 

• P = a general, specified point. 


2.4 COORDINATE SYSTEMS 

The following abbreviations will be used for the indicated coordinate 
systems. These coordinate systems are discussed in Applicable Document 1 
listed in Section 3.0. For all of them, the Y-axis forms a right hand 
orthogonal triad with the X and Z axes. 

• M50 - The Mean of 1950 is the Earth-centered inertial system whose X-axis 

is in the direction of the equinox at the beginning of the Bessel ian 
year 1950, and the X-Y plane is in the mean equator of that epoch. 

The Z-axis lies along the mean Earth's rotation axis of that epoch. 

• TOD - The True Of Date is the Earth-centered inertial system whose X-axis 

is in the direction of the Vernal Equinox of the midnight prior to 
launch. The XY plane is in the equitorial plane of that epoch, and 
the Z-axis lies along the Earth's rotation axis of that epoch. 

• ECI - The Earth-Centered Inertial system has its X-axis in the mean equator 

of the epoch (midnight prior to launch) and remains (i.e., is 
inertial) in the direction of the Greenwich meridian at that time. 

The Z-axis lies along the Earth's rotation axis of that epoch. 

t GEO - The GEOgraphic frame is Earth-centered and Earth-fixed with the X- 
axis lying in the equator of date and passing through the Greenwich 
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Note that h^ » Jl-? mutually orthogonal. By linear 

transformation of coordinates, each h^^ can be given in M50 by according 
to 



= (-simjj) + (cos<p) , and 
= (cos\|) cose) + (sinij) cos e) (-sine) • 

In accordance with foregoing arguments, the topodetic Euler angle rates 
are given by the projections of the body angular velocity _U as seen in the 
topodetic frame onto the axes » JI2 » il3 Euler angle 

rotations ip , e , and (p : 

i = U * Hi » 

® “ H * H2 ’ 

i = ]J • _H 3 . 


4. REFERENCE 

OFT Ascent/Descent Ancillary Data Requirements Document , Mission Planning and 
Analysis Division, Johnson Space Center, NASA, 78-FM-40, Rev. 2, JSC-14370, 
February 1982. 



meridian. Its Z-axis lies along the Earth's rotation axis of that 
date. 

• TOP - The TOPodetic system is an inertial system whose origin is on the 

Navigation Base with the X-axis pointing northward along the local 
meridian. The Z-axis points downward in the direction normal to the 
reference ellipsoid (presently, the Fischer Ellipsoid of 1960). Note 
that the topodetic frame is non-rotating ; a new frame is redefined at 
each time point. 

• BOD - The vehicle BODy system has its origin centered at the Navigation 

Base (not at the Center of Mass). The X-axis points toward the 
Orbiter's nose, and the Z-axis points down through the Orbiter's 
bottom. 

• PLM - The PLuMbline coordinate system is Earth-centered and inertial. Its 

axes are parallel to a right handed orthogonal system centered at the 
launch site and fixed inertially at the time of SRB ignition (i.e., 

/N 

GET=0). Define n as the unit vector normal to the tangent plane at 
the launch site; its direction is along the gravity gradient at that 
point. Also define a unit vector k at the launch site, in the 
tangent plane, and pointing in the direction of the launch azimuth. 

/V /s 

The unit vectors n and k are normal to each other and to the third 
unit vector, j , of the triad. Then the plumbline X, Y, and Z axes 
are parallel and positive in the direction of the unit 
vectors n, j, and k at the instant of SRB ignition. The Plumbline 
system is used only during ascent. 

• LF - The Landing Field (or Runway) coordinate system is a rotating 

coordinate system whose origin is at the intersection of the runway 
centerline and threshhold. The X-axis lies along the centerline, and 
is positive in the direction of rollout. The Z-axis is normal to the 
reference (i.e., Fischer) ellipsoid at the origin and is positive 
toward the Earth's center. 
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4.0 THE OUTPUT PRODUCTS 


The output product of the Ascent/Descent analysis program consists of two 
files, BETDATA and NAVBLK. Each is produced by a major program. The purpose 
of this section is to describe the methods by which these programs calculate 
the output product files. 


4.1 THE BETDATA FILE 

The BETDATA file consists of 239 double precision words calculated by the 
program OPIP and subroutines called by it. The following Section 4.1.1 
describes the processing within OPIP, and Section 4.1.2 describes the 
calculations of subroutines called by OPIP. 


4.1.1 OPIP Calculations 


The contents of the BETDATA file are discussed on the following pages in the 
order in which the words are calculated. The OPIP program and its subroutines 
calculate terms, each of which is represented by a name and an algebraic 
symbol and may consist of more than one double precision word in the file. 

For example, the M50 position vector has the term name M50, is designated by 
the algebraic symbol fg[M50], and consists of three double precision words. 

NB 

The following discussion provides the term name, algebraic symbol, the 
program, or subroutine in which the terms is computed, and the algorithm by 
which the term is computed. 
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BET OUTPUT PRODUCTS 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Time from SRB ignition (in seconds); i.e., Ground Elapsed Time (GET) 

1 GET OPIP Read next output time, from TRJATTDATA 

file 

OPIP Read SRB ignition time, from OPIP. IN 

file. 

OPIP tg^j = tpjgxT - ^SRBI 


Ground time in day : hour :min: sec of UT since the start of the year, but where day is 
the day number. That is, ground time is the ellapsed Astronomical UT since the 
start of the year + 1 day. 


2 

3 

4 

5 


GMTDAY 
GMTHR 
GMTMIN 
GMT SEC 


^NEXT 


TELAPSE Time is maintained internally in seconds 
from the astronomical UT start of year. 

The program is initialized with a Launch 
time (year, day number, seconds from the 
the UT midnight day of launch), stepped 
forward in time, and updated with time-tagged 
observations. Data records are written 
with GMTSEC as double precision. 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Onboard (Shuttle) time in day : hour :min: sec since SRB ignition 


6 SGMTDAY \ 

7 SGMTHR I 

8 SGMTMIN I 

9 SGMTSEC ’ 


OPIP 


Compute onboard time of next 
output, 


^NEXT “ ^t ^^drift * ^os 


'sNext 


1 - At 


drift 


= Difference (or error) between 
onboard and ground times. 


At = Drift rate of onboard 
clock. 


tos = time, in seconds, from midnight 

prior to launch as measured by the 
onboard clock. A value of is 
assumed. 


TELAPSE The output is computed in the same way as 
are words 2, 3, 4, and 5. The times are 
onboard times, however, rather than ground 
times (eg., replaces 


Year of present (Ground) time, UT 

10 GMTYR TELAPSE The present UT year is calculated along 

with words 2, 3, 4, and 5. 
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BET OUTPUT PRODUCTS (continued) 

OEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Position of the Nav Base relative to the (inertial) M50 coordinate frame 


11 

M50(l) 



Extracted from the state vector: 

12 

(2) 

i f, [M50] 

OPIP 

Tq [M50] = M50STAT(1,2,3) read in 



NB 


NB 

13 

(3) ^ 

f 


from the input file TRJATTDATA. 


Velocity of the Nav Base relative to the (inertial) M50 coordinate frame 


14 

DM50(1) s 


Extracted from the state vector: 

15 

(2) V r [M50] 

OPIP 

• 

[M50] = M50STAT(4,5,6) read 


1 NB 


NB 

16 

(3) ' 


in from the input file TRJATTDATA 


Contact (or sensed) acceleration of the Nav Base relative to the (inertial) M50 
coordinate frame 


17 

DDM50(1) X 


Extracted from the state vector: 

18 

(2) > r. [M50] 

OPIP 

Tq [M50] = M50STAT(7,8,9) read 


j NB 


NB 

19 

(3) ' 


in from the input, file TRJATTDATA 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Position of the Orbiter Center of Mass (CM) relative to the M50 coordinate frame 

20 M50C(1) . OP IP Read in r,^gCB0D] from TRJATTDATA file. 

I CM 

21 (2) \ r^ [M50] M3TX1 ^ ^NB 

I CM CM CM 

22 (3) ' VADD Tg [M50] = r^ [M50] + r,^g[M50] 

CM NB CM 


Velocity of the Orbiter CM relative to the M50 coordinate frame 

OPIP Read in illiCBOD] = M50STAT(13,14,15) 

Tg [M50] CROSS ^ 

CM CM CM 

M3TX1 ^ 

CM CM 

• • 

VADD Tg [M50] = Tg [M50] + v^g[M50] 

CM NB CM 


23 DM50C(1) 

24 (2) 

25 (3) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Sensed Acceleration of the Orbiter Center of Mass relative to the M50 coordinate 
frame 


26 

DDM50C(1) . 


OPIP 

Read in = 

= M50STAT(16,17,18) 

27 

(2) 

> r [M50] 

CROSS 

^.bCbod] = 

i!i[60D] X r^gCBOD] 



CM 


CM 

CM 

28 

(3) ^ 

1 

CROSS 

af^gCBOD] = 

i^EBOD] X v^g 





CM 

CM 

• 




VADD 


+ S[B0D] X P^gEBOD] 






CM 




M3TX1 

aNB^MSO] = 

[M50 ->■ BOD]^ 


CM CM 

VADD rg [M50] = r^ [M50] + a^g[M50] 
CM NB CM 


Position of 
29 GEO(l) 


30 (2) 

31 (3) 


the Nav Base in the GEOgraphic coordinate system 

M3X1 fg [ECI] = [M50 ^ ECI] rg[M50] 
NB NB 


Tg [GEO] 
NB 


ECIGEO 


[ECI ■*’ GE 03 ~ fctn(u)g, 


Pg [GEO] = [ECI ^ GEO] TgCECI] 
NB NB 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Velocity of the Nav Base in the GEOgraphics coordinate system 


32 DGEO(l) ^ 

I 

33 (2) > Tg [GEO] 

I NB 

34 (3) ' 


M3X1 


ECIGEO 


• • 

Pg [EC I] = [M50 ^ ECI] fg [M50] 

NB NB 

[ECI ->■ GEO] = fctn(o)g, 

uig = Earth's rotation rate vector 

• • 

Pg [GEO] = [ECI > GEO] [ECI] 

NB NB 

-itig X rg [GEO] 

NB 


Sensed Acceleration of the Nav Base in 

35 DDGEO(l) . M3X1 

36 (2) > aj^g[GE0] M3X1 

I sensed 

37 (3) ' 


the GEOgraphic coordinate system 

fg [ECI] = [M50 ^ ECI] fg [M50] 
NB NB 

a^g[GE0] = [ECI > GEO] r^ [ECI] 
sensed NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Total Acceleration of the Nav Base expressed in the GEOgraphic Coordinate System 


38 DDTGEO(l) 


39 


( 2 ) 


EARP 


r<g [GEO] 
NB 


Calculate the gravitation acceleration at 
the Nav Base in ECI frame. 


^NB 


[ECI] = fctn(r^ [ECI], 
NB 


40 


( 3 ) 


VADD 


K [ECI] = + [M50 


NB 


ECI] rg [M50] 
NB 


ECIGEO Tg [GEO] = [ECI > GEO] r^ [ECI] 

NB NB 

- X (dig X Tg [ECI]) 

NB 

- 2dig X [ECI ^ GEO] Tg [ECI] 

NB 


Gravitational Acceleration of the Nav Base expressed in the GEOgraphic coordinate 
system 

41 DDGGEO(l). EARP %CECI] = fctn(rg [ECI], t^^^^) 

I NB 

42 (2) \ Q^bCGEO] M3X1 = CECI ^ GEO] %[ECI] 

43 (3) ' 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Initial computations for all wind-relative velocities which follow 

SELMETA Read the meteorological data 
(Ascent) 

input file CMETDATA, interpolate 

DESCMET to the location ($q£qq. ^gEOD^ 
(Descent) and generate the array METDATA for that 
location. Define the following terms 
from METDATA: 


«s(w„) 

«n("h) 


= Horizontal wind speed 
= Systematic uncertainty in w^j 
= Noise uncertainty in w^ 




= Direction of horizontal wind 
(= 0° North and is positive 
clockwise from that) 


= Systematic uncertainty 


= Noise uncertainty in 

<^|y|(Wy) = Noise uncertainty in 

vertical wind speed 


The total uncertainty in horizontal wind 
speed, a^(w^), is 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


2 2 1/2 
Ol-(w^) = ] 

The total uncertainty in horizontal 

wind direction, is 




'N' WH' 


Wind-relative velocity of the Nav Base projected onto the GEOgraphic coordinate 


system 

44 DGEOW(l) X 


OPIP 

Set conversion constant Cq/j^ =degrees/radian 

45 (2) 

> V [GEO] 
wind 
NB 


Read R^^ Earth's equitorial radius and 
Rgp Earth's polar radius 

46 (3) ^ 



from input file OPIP. IN 


[GEO] is computed above in words (29-31) 
NB 


GEOGEOD Compute Geodetic latitude (I'q^qq)* 
longitude (-^q£ 0 D^’ height (hg^QD^’ 
and declination, 6 . (See description of 
GEOGEOD in Section 4.1.2) 

GEOTOP [GEO ^ TOP] = i^ctn(Ag^Qp, 4 .^^^^) 

• • 

M3X1 Tg [TOP] = [GEO > TOP] [GEO] 

NB NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DGEOW continued 

WINDTOP ~ 

$(^wi ~ fctn(0(W|^), w^, o(W|_j)) 

(See the description of WINDTOP in Section 
4.1.2 for the expressions for the matrix 
elements in (|:(v^^-^j[T0P]) . ) 


vsub 

Q- 
1 — 1 
—1 

0 

1 1 

It 



NB 

NB 

M3TX1 

"wlndCSEO: - 

[GEO » TOpf 5„,„<([T0P] 


NB 

NB 


Position of the Center of Mass in the GEOgraphic coordinate system 


47 

GEOC(l) V 

OPIP 

Read 

[M50 ->» ECI] from OPIP. IN 

48 

(2) > r. [GEO] 

OPIP 

Read 

r.g[B0D] and [M50 ^ BOD] 


1 CM 



CM 

49 

(3) ' 


from 

TRJATTDATA 


ECIGEO fg [GEO] defined previously in 
NB words (29-31) 


M3TX1 ^ 

CM CM 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

GEOC continued 


M3X1 


ECIGEO 


VADD 


ri^gCECI] = [M50 ^ ECI] r^g[M50] 

CM CM 

r^gCGEO] = [ECI ^ GEO] Tf^gCECI] 

CM CM 

[GEO] is defined in words (29-31) 
NB 

[GEO] = Tg [GEO] + r^g[GEO] 

CM NB CM 


Velocity of the Center of Mass in the GEOgraphic coordinate system 


50 

DGEOC(l) 

OPIP 

Read in w[BOD] = M50STAT(13 ,14,15) 

51 

1 • 

(2) y r. [GEO] 

CROSS 

• 

r^g[B0D] = i!;[BOD] X r^g[B0D] 


1 CM 


CM CM 

52 

(3) ^ 

M3TX1 

• • 

rfjg[M50] = [M50 ^ BOD]^ 




CM CM 



M3X1 

• • 

r,^g[ECI] = [M50 > ECI] r^g[M50] 


CM CM 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DGEOC continued 


ECIGEO 


VADD 


r^gf^^EO] = [ECI ^ GEO] Pf^gCECI] 
CM CM 



Sensed acceleration of the Center of Mass in the GEOgraphic coordinate system. 

OPIP Read in ^^>[800] = M50STAT(16,17,18) 


a^^n^CGEO] CROSS " ^[BOD] X r^g[BOD] 

sensed CM CM 

CROSS Centrifugal acceleration: 

acENCBOD] = i[BOD] X r^gCBOD] 

CM 

CROSS Rotational acceleration: 

agoT-CBOD] = ;:5[B0D] X r^g[BOD] 

CM 

VADD ^ROT>^®®'^^ 

CM 

MT3X1 ^ 

CM CM 


53 DDGEOC(l) 

54 (2) 

55 (3) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DDGEOC continued 


M3X1 


ECIGEO 


rNR[ECI] = [M50 > ECI] rj^gCMSO] 

CM 


NB‘ 

CM 


ac„[GEO] 

sensed 


[ECI > GEO] 

CM 

CM 

- X [ECI ^ GEO] rf^gCECI] 

CM 


Wind-relative velocity of the Center of Mass projected into the Geographic system 


56 DGEOWC(l) 


OPIP v[GEO] computed previously, words (44-46) 
wind 
NB 


57 


58 


( 2 ) 


( 3 ) 


v[GEO] 

wind 

CM 


rj(|g[GE0] calculated previously, in order 
CM 

to calculate r^ [GEO], words (50-61) 

CM 


VADD 


V [GEO] 
wind 
CM 


V [GEO] + r,^g[GE0] 
wind 1,^ 

NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Position of the Nav Base relative to the landing runway expressed in Landing Field 
coordinates 

59 LF(1) . 6E0LF 


60 (2) \ \ [LF] 

1 NB 

61 (3) ^ OPIP fg [GEO] computed previously, words (29-31) 

NB 

VSUB \ [GEO] = r^ [GEO] - Pg[GE0] 

NB NB L 

M3X1 r^ [LF] = [GEO > LF] f|_ [GEO] 

NB NB 


First Pass obtain [GEO >■ LF] and 
?g[GE0] from Landing site geophysical 
L 

parameters 


Velocity of the Nav Base relative to the landing runway in Landing Field coordinates 


62 DLF(l). 

63 (2) \ v^ [LF] 

1 NB 

64 (3) ' 


OPIP 


M3X1 


iTg [GEO] computed previously, words (32-34) 
NB 

V, [LF] = [GEO > LF] r [GEO] 

NB NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Sensed acceleration of the Nav Base relative to landing runway in Landing Field 


coordinates 




65 DDLF(l) X 


OPIP 

Tg [GEO] computed previously, words (38-40) 
NB 

66 (2) j 

[LF] 

i NB 

M3X1 

a, [LF] = [GEO ^ LF] f. [GEO] 
NB NB 

67 (3) ^ 





Geodetic velocity of the Nav Base projected in the TOPodetic coordinate frame 


68 DTOP(l) K 

69 (2) V Vg [TOP] 

1 NB 

70 (3) ' 


GEOGEOD Geodetic latitude and longitude, 

ij) and X , respectively, were 
GEOD GEOD 

computed previously (words 44-48) from 
the geophysical parameters and 
fg [GEO] which was computed 
NB 

previously, words 29-31 


GEOTOP [GEO ^ TOP] = fctn(Xg^Qp, 

M3X1 Vg [TOP] = [GEO > TOP] [GEO] 
NB NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Geodetic sensed acceleration projected in the TOPodetic coordinate frame 


71 

DDTOP(l) X 


M3X1 

[TOP] = [GEO > TOP] 

72 

J 

> ag [TOP] 


NB 

sensed 

73 

(3) ^ 

NB 

' sensed 




rg [GEO] 
NB 


Uncertainty (one sigma) in Vg [TOP] 

NB 

74 DTOPU(l) K OP IP 

75 (2) > o(Vg [TOP]) OPIP 

I nb 

76 (3) ' M3X3UN2 

M3X3UN2 

M3X3UN2 

SIGMAS 


Read (|:(§[M50]) from input file TRJATTDATA 


f(rg [M50]) is extracted from (t:($[M50]) 

NB 

• • 

(f(fg[ECI]) = [M50 > ECI] <|:(fg[M50]) [M50 > ECI]^ 
NB NB 

• • 

<P(rg[GEO]) = [ECI ^ GEO] (t:(fQ[ECI]) [ECI - GEof 

NB NB 

$(Vg[T0P]) = [GEO > TOP] (|:(fg[GE0]) [GEO TOP]^ 
NB NB 

5(»V™P]) = hi'"’' 

NB 

where [C.j] = (I:(Vg[T0P]) 

NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 

Uncertainty (one sigma) in a„ [TOP] 


NB 

sensed 


77 DDTOPU(l) 


OPIP (|:(r[M50] is extracted from <1;($[M50]) 


78 


79 


( 2 ) 


( 3 ) 


5(a„ [TOP]) M3X3UN2 


NB 

sensed 


ll:(r[ECI]) = [M50 > ECI] ([(fCMSO]) [M50 > ECI] 


M3X3UN2 (p(r[GEO]) = [ECI ^ GEO] (t(?[ECI]) [ECI > GEO]^ 


M3X3UN2 {f(a[TOP]) = [GEO > TOP] (|:(r[GEO]) [GEO > TOP]^ 


1/2 1/2 i /2 T 

SIGMAS a(a^ [TOP]) = [C^^ ’^22’ ^33 ^ 

NB 

sensed 

where [C.j] = (f(ag[TOP]) 

NB 

sensed 


Wi nd-rel ati ve 

80 DTOPW(l) 

81 (2) 

82 (3) 


velocity of 


V [TOP] 
wi nd 
NB 


the Nav Base projected on to the TOPodetic frame 


WINDTOP V [TOP] = [-w^ cos 6^^, -w^ sin of 


WINDTOP $(v [TOP]) = fctn(a(wj^), cj(0|^^), w^. o(wJ) 
NB 

(See description of WINDTOP in Section 4.1.2 
for the expressions for the matrix elements) 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Velocity of the Nav Base in the GEOgraphic coordinate system, but projected into the 
vehicle BODy axis coordinate system 

ECIGEO Tg [GEO] computed previously, words (32-34) 
NB 

Vg[B0D] M3X3T [ECI -»■ BOD] = [M50 ^ BOD] [M50 ♦ ECI]^ 

NB 

M3X3T [GEO + BOD] = [ECI ^ BOD] [ECI > GEO]^ 

M3X1 Vg [BOD] = [GEO ^ BOD] r^ [GEO] 

NB NB 


83 DBOD(l) 

84 (2) 

85 (3) 


Sensed acceleration of the Nav Base projected onto the vehicle BODy system 


86 

DDBOD(l) ^ 

OPIP 

Read r^ [M50] and [M50 BOD] 


1 


NB 


I 


from the input file TRJATTDATA 

87 

(2) > L [BOD] 




1 NB 


.. 


1 

M3X1 

[BOD] = [M50 ^ BOD] r^ [M50] 

88 

(3) ' 


NB NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Wind vector projected into the vehicle BODy coordinate system 


89 WRODV(l) 


WINDTOP Vg [TOP] = C-Wy cos0,,^J, sine,^,^, 0]‘ 

wind 


90 


(2) > Vq [BOD] M3X1 Vg [BOD] = [TOP ^ BOD] v^ [TOP] 
wind wind wind 


91 


( 3 ) 


Wind-relative velocity of the Nav Base expressed in vehicle BODy coordinates 


92 DBODW(l) 


OPIP Vg [TOP] computed previously, words (68-70) 

NB 


93 


94 


( 2 ) 


( 3 ) 


V [BOD] 
wind 
NB 


WINDTOP Vg [TOP] = [-w^ "”h 

wind 


VSUB 


M3X1 


V [TOP] 
wind 

NB 

V [BOD] 
wind 

NB 


Vg [TOP] - Vg [TOP] 
NB wind 

[TOP ^ BOD] V [TOP] 
wi nd 
NB 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Uncertainty of velocity of Nav Base in GEOgraphics coordinates projected into BOU 
system 


95 DB0DU(1)\ 


OPIP (|:(rg [M50]) extracted from (f(§[M50]) 

NB 


96 


(2) I a(Vg [BOD]) OPIP 
I NB 


>1^) extracted from $(§[M50]); 
(|) = roll, 9 = pitch, = yaw 


97 

I 

I 

i 


( 3 ) 


M3X3UN2 (f(rg [ECI]) = 

NB 

[M50 -»■ ECI] (|:(rg [M50]) [M50 > ECI]"'^ 
NB 

ECIGEO [ECI ■>■ GEO] computed previously, for 
words (19-31) 

M3X3UN2 (p(rg [GEO]) = 

NB 

[ECI ^ GEO] [ECI]) [ECI > GEO]^ 
NB 

M3X3T [ECI ^ BOD] = [M50 ^ BOD] [M50 > ECI]^ 


M3X3T [GEO ^ BOD] = [ECI ^ BOD] [ECI ^ GEO]^ 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DBODU continued 

M3XUN2 $'(Vg [BOD]) = 

NB 

[GEO ^ BOD] (pCrg [GEO]) [GEO > BOD]^ 
NB 


SBCALC 


(1:(Vq [BOD]) 
NB 



(See the description of SBCALC in 
Section 4.1.2) 


SIGMAS S(vg [BOD]) = 

NB 

where [C.^] = ^(Vg [BOD]) 
NB 


Uncertainty of sensed acceleration of Nav Base projected onto the BODy system 


98 

DDBODU(l) ^ 

OPIP 

<f(r| [M50]) & (p( 4 ), 0, 4 <) from (f(![M50]) 



1 

NB 

99 

(2) 

[BOD]) OPIP 

[M50 ^ BOD] from TRJATTDATA 



NB 


100 

1 

(3) > 

M3X3UN2 

[BOD]) = 


NB 

[M50 ^ BOD] (|:(rQ [M50]) [M50 > BOD]"^ 
NB 
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BET OUTPUT PROeUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DDBODU continued 


SPBCALC 


$(ag [BOD]) = [BOD]) + ( 

NB NB 


effects of body 
axis errors 


1 


(See the description of SPBCALC in 
in Section 4.1.2) 


SIGMAS $(ag [BOD]) = [Cl//^ 0^3 C33^/^]'^ 
NB 

where [C.^] = (|;(ag [BOD]) 

NB 


Uncertainty in wind 

vector projected 

onto BODy axis 

coordinates 

101 

WBODVU(l) ^ 

OPIP 

$(<!’. 9 ) extracted from (f(^[M50]) 

102 

(2) > a(; 

/g [BOD]) OPIP 

[M50 ^ BOD] read from the TRJATTDATA file 



wind 



103 

(3)^ 

OPIP 

[M50 > ECI] 

read from the OPIP. IN file 



M3X3T 

[ECI > BOD] 

= [M50 > BOD] [M50 > ECI]^ 



GEOTOP 

[GEO -»■ TOP] 

fctn('l>QEOD’ ^GEOD^ 



ECIGEO 

[ECI -»• GEO] 

= fctn(o)g, t, Og [ECI], 





NB 

• • 





rg [ECI], Tg [ECI], Tg [ECI]) 





NB NB NB 



M3X3 

[ECI > TOP] 

= [GEO ^ TOP] [ECI ^ GEO] 



M3X3T 

[TOP ^ BOD] 

= [ECI ^ BOD] [ECI ^ TOP]^ 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

WBODVU continued 

WINDTOP (f{v [TOP]) = fctn(a(Wj^), a(e^^^), w^, a(w^)) 
wi nd 

(See description of WINDTOP in Section 4.1.2 
for expressions for the matrix elements) 

M3X3UN2 (f'(v [BOD]) = 
wi nd 

[TOP -»■ BOD] (j:(v [TOP]) [TOP > BOD]^ 
wi nd 


SPBCALC (p(v [BOD]) = f (v [BOD]) + {aIIrerrors‘'°'^^ 

dAibtrrrurs 

wind wind 

(See description of SPBCALC in Section 
4.1.2) 

SIGMAS a(Vg [BOD]) = ' ^33^^^^^ 

wind 

where [C.^] = $(Vg [BOD]) 
wind 


Uncertainty in the wind-relative velocity of the Nav Base expressed in BODy 
coordi nates 


104 


105 


106 


DBODWU(l) 


( 2 ) 

( 3 ) 


?(v [BOD]) 
wind 
NB 


OPIP 


OPIP 


OPIP 


$(<t'» 9» 'p) extracted from (f(§[M50]) 

[TOP -»■ BOD] computed previously, 
words (101-103) 

(p(Vg [TOP]) computed previously, 

NB 

words (74-76) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DBODWU continued 


WINDTOP V [TOP] = [-w^ cos 6^^, -w^ sin 0]^ 

wind 

WINDTOP (Kv [TOP]) = fctn(o(w^), o(b^^), w^^, u(w^)) 
wi nd 

(See description of WINDTOP in Section 
4.1.2 for the expressions for the matrix 
elements) 


VSUB 


OPIP 


M3X3UN2 


SPBCALC 


V [TOP] = Vg [TOP] - V [TOP] 
wind wind 

NB 


[TOP]) = (l:(Vg [TOP]) + Uw [TOP]) 
wind .,o wind 

NR 


(p(v [BOD]) = 
wi nd 
NB 

[TOP ^ BOD] (p(v [TOP]) [TOP > BOD]' 


wi nd 
NB 


Uy [BOD]) = (p(v [BOD]) + { 


wind wind 

NB NB 


effects of body 
axis errors 


1 


(See description of SPBCALC in Section 
4.1.2) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 


word name symbol 

DBODWU continued 


program algorithm 


SIGMAS 


+ 1/2 1/2 1/2,T 

a(v [BOD]) = [C,, > ^zz > ^33 ^ ^ 

wind 
NB 

where [C. •] = (|:(v [BOD]) 
wi nd 
NB 


Quaternion for transforming from M50 coordinates of BODy coordinates 


107 050(1) 

108 (2) 

109 (3) 

no (4) 


OPIP 

q [M50 ^ BOD] 


Read directly from the input file TRJATTDATA 
and equivalenced to the output buffer array 
elements 


Index of selected IMU 

111 PIMU OPIP Read directly from the input file TRJATTDATA 

and equivalenced to the output buffer array 
element 


Transformation from the ECI to launch-site-located Plumbline coordinate system 

ECIPLM [ECI > PLM] is calculated on the first 
pass for ascent (only) from geophysical 
parameters, launch site parameters, and 
the time of SRB ignition 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Position of the Nav Base relative to the Earth's center expressed in Plumbline 
coordinates (Ascent only) 


112 PLM(l) 


OPIP Read [M50 > ECI] from input file OPIP.IN 


113 (2) 


fg [PLM] OPIP 
NB 


Read [M50] from the input file TRJATTDATA 


NB 


114 (3) 


M3X1 r. [ECI] = [M50 ECI] r. [M50] 


NB 


NB 


M3X1 r. [PLM] = [ECI > PLM] [ECI] 


NB 


NB 


Velocity of the Nav Base relative to the Earth's center expressed in Plumbline 
coordinates (Ascent only) 


115 

DPLM(l) . 

OPIP 

• 

Read L [M50] from TRJATTDATA 




NB 


116 

1 

(2) V r^ [PLM] 

M3X1 

• 

rg [ECI] = [M50 ^ ECI] 

• 

Tg [M50] 


1 NB 


NB 

NB 

117 

(3) ^ 

M3X1 

• 

Tg [PLM] = [ECI ^ PLM] 

• 

Pg [ECI] 




NB 

NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Sensed acceleration of the Nav Base relative to Earth's center expressed in 
Plumbline coordinates (Ascent only) 


118 

r)DPLM(l) \ 


OPIP 

Read r^ 

[M50] from TRJATTDATA 





NB 


119 

(2) 

1 

f, [PLM] 


1 1 

o 

UJ 
( 1 

: 

= [M50 ^ ECI] Tq [M50] 


1 

NB 


NB 

NB 

120 

(3) ) 



[PLM] 

= [ECI > PLM] fg [ECI] 





NB 

NB 


Uncertainty in the position of the Nav Base relative to the Earth's center 
expressed in Plumbline coordinates (Ascent only) 


121 PLMU(l) 

122 ( 2 ) 

123 (3) 


o{r^ [PLM]) 
NB 


OPIP 


M3X3UN2 


(f(rA [M50]) is extracted from (|:(§[M50]) 


NB 


<t:(rg [ECI]) = 

NB 

[M50 > ECI] l|:(rg[M50]) [M50 
NB 


ECI]' 


M3X3UN2 [PLM]) = 

NB 

[ECI ^ PLM] [ECI]) [ECI > PLM]"^ 
NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

PLMU continued 


SIGMAS 


o(rg [PLM]) = C 22 ^/^ 033 '/^]^ 


NB 


where [C.^] = Ur ^ [PLM]) 
NB 


Uncertainty in the velocity of the Nav Base relative to the Earth's center 
expressed in Plumbline coordinates (Ascent only) 


124 OPLMU(l) . 

I 

125 ( 2 ) V o(rg [PLM]) 

I NB 

126 (3) ^ 


OPIP 


M3X3UN2 


M3X3UN2 


SIGMAS 


U^Q [M50]) is extracted from (])(^[M50]) 

NB 

|(rg [ECI]) = 

NB 

[M50 > ECI] Ur^ [M50]) [M50 ^ EC if 
NB 

(f(r0[PLM]) = [ECI > PLM] (l:(rg[ECI]) [ECI > PLMf 
NB NB 

o(fg, [PLM]) = ^ 33 '/^]^ 

NB 

where [C.^] = (Krg[PLM]) 

NB 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Uncertainty in the sensed acceleration of the Nav Base relative to the Earth's 
center expressed in Plumbline coordinates (Ascent only) 

<^(r^ [M50]) is extracted from (|;(^[M50]) 

NB 

(Krg [ECI]) = 

NB 

[M50 ^ ECI] ({;(rg [M50]) [M50 ^ ECI]^ 
NB 


127 DDPLMU(l) 


OPIP 


128 (2) 


o(rg [PLM]) M3X3UN2 
NB 


129 (3) 


M3X2UN2 (l:(rg [PLM]) = 

NB 

[ECI ^ PLM] [ECI]) [ECI ^ PLM]^ 
NB 


SIGMAS o(rg [PLM]) = [C^^ > ^22 ' ^33 ^ 

NB 

where [C.^] = [PLM]) 

NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Wind-relative velocity of the Nav Base but projected onto Plumbline coordinates 
(Ascent only) 


130 

PPLMW(l) \ 


VSUB 

131 

(2) 

V V [PLM] 

M3X3T 


wi nd 




NB 

ECIPLM 

132 

(3) y 




[TOP] was computed previously for 
wi nd 
NB 

words (104-106) 

[ECI + TOP] was computed previously for 
words (101-103) 

[ECI + PLM] was computed previously just 
before words (112-114) 


M3TX1 


M3X1 


V [ECI] 
wi nd 
NB 


V [PLM] 
wind 
NB 


[ECI ^ TOP]^ V [TOP] 
wi nd 
NB 


[ECI + PLM] V [ECI] 
wind 
NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Uncertainty in the wind-relative velocity of the Nav Base projected onto the 
Plumbline coordinate system (Ascent only) 


133 


134 



WINDTOP 


M3X3UN2 


(|;(v [TOP]) computed previously for 
wind 

words (80-82) 

$(Vg [TOP]) computed previously for 
NB 

words (74-76) 


135 (3) 


OPIP (l;(v [TOP]) = (Kv [TOP]) + (p(v [TOP]) 
wi nd wi nd 

NB 


M3X3UN1 (j:(v [ECI]) = 
wi nd 
NB 

[ECI > TOP]' $(v [TOP]) [ECI > TOP] 
wind 
NB 


M3X3UN2 (I:(v [PLM]) = 

wind 
NB 

[ECI > PLM] [ECI]) [ECI > PLM]' 
wind 
NB 

SIGMAS S(v [PLM]) = [C,//^ C„^/^ 

wind 
NB 

where [C. .] = (|:(v [PLM]) 

^ wind 

NB 


40 



BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Position of the vehicle Center of Mass relative to Earth's center expressed in 
Plumbline coordinates (Ascent only) 



\ [PLM] 
CM 


OPIP 

ECIPLM 


M3X1 


Read f|^g[B0D] from the input file TRJATTDATA 
CM 

[ECI PLM] computed previously before 
words (112-114) 

fj^eCMSO] = [M50 ^ BODf r^gCBOD] 

CM CM 


M3X1 


M3X1 


f^gCECI] = [M50 > ECI] rfjg[M50] 
CM CM 

r^gCPLM] = [ECI ^ PLM] r^^gEECI] 
CM CM 


OPIP fg [PLM] is word set (112-114) 
NB 

VADD fg [PLM] = Tg [PLM] + r,^g[PLM] 
CM NB CM 



BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Velocity of the vehicle Center of Mass relative to the Earth's center expressed in 
Plumbline coordinates (Ascent only) 


139 DPLMC(l) 


CROSS v.qCBOD] = i[B0D] X CmdCBOD] 


CM 


NB‘ 

CM 


140 


(2) > fg [PLM] M3TX1 = [M50 > BOof v^g[B0D] 


CM 


CM 


CM 


141 (3) 




CM 


CM 


M3X1 = [ECI ^ PLM] rf^gCECI] 


CM 


CM 


OPIP r» [PLM] is word set (115-117) 


NB 


VADD Tj) [PLM] = Tg [PLM] + r^g[PLM] 


CM 


NB 


CM 


Sensed acceleration of the vehicle Center of Mass expressed in Plumbline coordinates 
(Ascent only) 


142 

DDPLMC(l) 

143 

(2) 

144 

(3) 


Cg [PLM] 
CM 


CROSS 

VADD 

M3TX1 


a^g[B0D] = u[B0D] X v^g[B0D] 

CM CM 

+ w[B0D] X r,^g[B0D] 

CM 

r^g[M50] = [M50 > BOD]''' a^g[B0D] 
CM CM 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DDPLMC continued 


M3X1 = [M50 ^ ECI] ?j^g[M50] 

CM CM 

M3X1 = CECI ^ PLM] f,^gCECI] 

CM CM 

OPIP Tg [PLM] is word set (118-120) 
NB 

VADD Pg [PLM] = r^ [PLM] + 

CM NB CM 


Wind-relative velocity of the vehicle Center of Mass expressed in Plumbline 
coordinates (Ascent only) 

OPIP 


V [PLM] 
wind 
CM 

M3X1 


M3X3 


145 DPLMCW(l) 

146 (2) 

147 (3) 


V [TOP] computed previously for 
wind 
NB 

words (80-82) 

rj^g[PLM] computed previously for 
CM 

words (139-141) 

[ECI > TOP] computed previously for 
words (101-103) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

DPLMCW continued 

ECIPLM [ECI >■ PLM] computed previously for 
words (112-114) 


M3TX1 V [ECI] = [ECI > TOP]"^ v [TOP] 
wind wind 

NB NB 


M3X1 V [PLM] = [ECI > PLM] v [ECI] 

wind wind 

NB NB 

VADD V [PLM] = V [PLM] + rMo[PLM] 

wind wind 

CM NB 


Yaw, and its uncertainty, of the vehicle body with respect to a wind-relative local 


frame 
148 YAW 

<1^ 

EULANG 

Given in the description of EULANG 

149 YAWU 



in Section 4.1.2 


Pitch, and its uncertainty, of the vehicle body with respect to a wind-relative 
local frame 


150 

PITCH 

6 

EULANG 

Given in the description of EULANG 
in Section 4.1.2 

151 

PITCHU 

a(e) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 


word 

name 

symbol 

program 

al gorithm 


Roll , 

and its 

uncertainty. 

of the vehicle body with respect to a 

wind-relative 

frame 






152 

ROLL 


EULANG 

Given in the description 

of EULANG 





in Section 4.1.2 


153 

ROLLU 

a(4>) 




Yaw Rate, and 

its uncertainty, of the 

vehicle body with respect to a wind-relat 

local 

frame 





154 

YAWD 


EULANG 

Given in the description 

of EULANG 





in Section 4.1.2 


155 

YAWDU 

a(i) 




Pitcin 

1 Rate, and its uncertainty, of the vehicle body with respect 

. to a 

wi nd- 

rel ati ve 

local frame 




156 

PITCHD 

0 

EULANG 

Given in the description 

of EULANG 





in Section 4.1.2 


157 

PITCHDU 

a(0) 




Roll 

Rate, and 

its uncertainty, of the 

vehicle body with respect 

to a 

wind- 

rel ati ve 

local frame 




158 

ROLLD 

• 

EULANG 

Given in the description 

of EULANG 





in Section 4.1.2 


159 

ROLLDU 

a(4i) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Bank angle 

160 BA $ BANK Given in the description of BANK 

in Section 4.1.2 


Longitude of the Nav Base 

161 GEOLONG ^GEOD GEOGEOD Given in the description of GEOGEOD 

in Section 4.1.2 


Geodetic Latitude of the Nav Base 

162 GEODLAT '*’gE 0D GEOGEOD Given in the description of GEOGEOD 

in Section 4.1.2 


Geodetic altitude (height) above the (Fischer) reference ellipsoid 
of the Nav Base 

163 GEODALT *^GE0D GEOGEOD Given in the description of GEOGEOD 

in Section 4.1.2 


Uncertainty in the geodetic altitude (height) of the Nav Base 

164 GEOHU ^^*^GE0D^ local vertical unit vector, n, which 

is normal to the reference ellipsoid be 

n = [n^, n^, n^] where 

^1 " ‘^°^%0D ^°^^GE0D 
^2 " ^°^‘*’gE0D ^^^^GEOD 
^3 "" ^^'^‘*’gE0D 


46 



BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

6E0HU continued 

OP IP Read M50C0V from input file TRJATTDATA 


OPIP 


M3X3UN2 


$(rg [M50]) = MVCOV(i.j); i = 1,3; j = 1,3 
NB 

$(rg [ECI]) = 

NB 

[M50 > ECI] (t;(rg [M50]) [M50 > ECI]^ 


NB 


M3X3UN2 $(rg [GEO]) = 

NB 

[ECI ^ GEO] (F(fg [ECI]) [ECI ^ GEO]^ 
NB 


OPIP 


®^^GEOD 


= n $(rg [GEO]) n"^ 
NB 


OPIP 


'^^^GEOD^ 


2 1 /2 


Geodetic altitude rate of the Nav Base 

165 DH *^GlEOD ^0 [GEO] computed previously, words (32-34) 

NB 

^GEOD^ ^ 

NB 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Uncertainty in geodetic altitude rate of the Nav Base 
166 DHU °(V,E0n^ 


<I:(r@ [M50]) = MVC0V(i+3J+3); i = L 
NB 


M3X3UN2 (l:(rg [ECI]) = 

NB 

[M50 ^ ECI] [M50]) [M50 
NB 

M3X3UN2 UK CGEO]) = 


NB 


[ECI ^ GEO] l(:(rg [M50]) [ECI 


NB 


o(hgE00^ = n Ur^ [GEO]) n^ 


NB 


‘’^^GEOD^ '^°('^GEOD^ ^ ^ 


Declination of the Nav Base 


167 DELTA 6^^ 


GEOGEOD Let r^ [GEO] = [x^. y^, ; then, 

NB 


6.,d = arctan ^ 

m 2 2 1/2 

(><5 - Vg ) ' 


3; j = 1,3 


ECI]^ 


GEO]^ 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Magnitude of the geocentric position vector of the Nav Base 

i 

168 RM Rg FPANG2 Rg = |fg [M50] | 

j NB NB NB 

I 

j 

Magnitude of inertial velocity of the Nav Base 

169 VM V^g FPANG2 V^g = 1?^ [M50]| 

NB 


Flight path angle as observed in the ECI coordinate frame 

170 FPM a FPANG2 Given in the description of FPANG2 

in Section 4.1.2 


Azimuth of the velocity observed in the ECI coordinate frame 

171 AZM A^ FPAN62 Given in the description of FPANG2 

of Section 4.1.2 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Slant range of the Nav Base measured from the launch site (Ascent) or Landing 


Site (Runway threshold, Descent) 

172 RS Rg OPIP 

NB 


GEODGEO 

VSUB 

OPIP 


Read geodetic position parameters 

(s). X (s), h (s)] 

GEOD GEOD GEOD 

for s = Origin of Launch Site (Ascent) 
or s = Origin of Landing Site (Descent), 
from input file OPIP. IN 

fjGEO] = fctn(4> (s), X (s), h (s)) 

® GEOD GEOD GEOD 

(Algorithm given in the description 
of GEODGEO in Section 4.1.2) 

r^ [GEO] = Tg [GEO] - rg[GE0] 

NB NB s 

Rj = IPj [GEO] I 
NB NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 

Transformation matrix, M50 frame to vehicle BODy frame 


173 

M502B0D(1,1) \ OP IP 

Read from the input file TRJATTDATA 

174 

(1.2) I 

and loaded directly into the output 

175 

(1.3) 1 

array. 

176 

(2.1) 1 


177 

(2,2) \ [M50 ^ BOD] 


178 

(2.3) ( 


179 

(3.1) 1 


180 

(3.2) 1 


181 

(3,3)/ 



Wind-relative velocity magnitude of the vehicle Center of Mass 

182 VTOP V . . FPAN6 Given in the description of FPANG 

wi nd 

CM in Section 4.1.2 


Uncertainty in the wind-relative velocity magnitude of the vehicle Center of Mass 

183 VTOPU ^^'^wind^ FPANG Given in the description of FPANG 

CM in Section 4.1.2 


Wind-relative flight path angle of the vehicle Center of Mass 

184 GAMTOP ^wind FPANG Given in the description of FPANG 

CM in Section 4.1.2 
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BET OUTPUT PRODUCTS (continued) 

DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Uncertainty in the wind-relative flight path angle of the vehicle Center of Mass 

185 6AMT0PU FPANG Given in the description of FPANG 

CM in Section 4.1.2 

Wind-relative azimuth of the vehicle Center of Mass 

186 PSITOP **^wind FPANG Given in the description of FPANG 

CM in Section 4.1.2 


Uncertainty in the wind-relative azimuth of the vehicle Center of Mass 

187 PSITOPU o(i|» . .) FPANG Given in the description of FPANG 

'wind' 

CM in Section 4.1.2 


Earth surface range from subvehicle (Nav Base) point to s, where s is the Launch 
Site origin (for Ascent analysis) or Runway threshold (for Descent analysis) 

188 S Rj^g^ OPIP Read geodetic position parameters, 

s [t (s), X (s), h (s)] 

GEOD GEOD GEOD 

for s = Origin of Launch Site (Ascent) 
or s = Origin of Landing Site (Descent), 
from input file OPIP. IN 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

S continued 


GEODGEO r^CGEO] = fctn(<)) (s), A (s), h (s)) 
® GEOD GEOD GEOD 

(Algorithm given in the description 
of GEODGEO in Section 4.1.2) 


Let r^CGEO] = [x^ , y^, 
s 

and r^ [GEO] = yj^g, Zj^g] 

NB 

Then 

2 2 2 2 
OPIP ^s ^s 

s 

^NB ^NB ^NB* 

NB 


OPIP 


OPIP 


‘^NBS ^ 


s s 

B is the angle, in radians, between the 
geocentric vectors to s and to the Nav Base 
and is calculated from 

2 2 2 


3 


= arc tan 


^0 ^0 " ^^0 ^0 ) 
s NB s NB 

R® R® cos 0 


where 


s NB 


Rg Rg cos 3 = TgCOEO] • Tg [GEO] 
s NB s NB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Wind-relative angle of attack for the vehicle Center of Mass 

189 ALPHA AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Uncertainty in the wind-relative angle of attack for the vehicle Center of Mass 

190 ALPHAU AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative guidance sideslip angle for the vehicle Center of Mass 

191 BETAP AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Uncertainty in the wind-relative guidance sideslip angle for the Center of Mass 

192 BETAPU ®(*^WG^ AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative dynamic pressure for the total Center of Mass velocity, in Ib/ft^ 

193 QAERO q AERODYN Given in the description of AERODYN 

in Section 4.1.2 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Uncertainty in the dynamic pressure for the velocity of the Center of Mass 

194 QAEROU a(q) AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative pitch dynamic pressure for the Center of Mass, in 1b-deg/ft^ 

195 QALPHA AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative yaw dynamic pressure for the Center of Mass, in Ib-deg/ft^ 

196 QBETA q^ AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative Mach number at the Center of Mass 

197 MACH M AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Uncertainty in the wind-relative Mach number at the Center of Mass 

198 MACHU a(M) AERODYN Given in the description of AERODYN 

in Section 4.1.2 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Wind-relative viscous parameter affecting the motion of the Center of Mass 

199 VOO V AERODYN Given in the description of AERODYN 

in Section 4.1.2 

Uncertainty in the wind-relative viscous parameter for the Center of Mass 

200 VOOU a(v) AERODYN Given in the description of AERODYN 

in Section 4.1.2 

Ambient atmospheric temperature in °R 

201 METDATA(3) T^ OPIP Placed in output buffer directly 

Uncertainty in ambient atmospheric temperature in ®R 

202 METDATA(7) Placed in output buffer directly 

Ambient atmospheric pressure in Ib/ft^ 

203 METDATA(4) P^ OPIP Placed in output buffer directly 
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BET OUTPUT PRODUCTS (continued) 

OEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Uncertainty in ambient atmospheric pressure in Ib/ft^ 

204 METDATA(8) Placed in output buffer directly 


O 

Atmospheric mass density in slugs/ft 

205 METDATA(5) OPIP Placed in output buffer directly 


Wind-relative equivalent air speed at the Center of Mass 

206 EAS AERODYN Given in the description of AERODYN 

CM in Section 4.1.2 


Uncertainty in the wind-relative equivalent air speed at the Center of Mass 

207 EASU AERODYN Given in the description of AERODYN 

CM in Section 4.1.2 


Load factor at the Center of Mass, in g's [ft/sec^] 

208 L n AERODYN Given in the description of AERODYN 

in Section 4.1.2 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Wind-relative drag over mass for the Center of Mass 

209 DOM *^D/M AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative lift over drag for the Center of Mass 

210 LOD Lq AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Wind-relative aerodynamic side slipangle for the Center of Mass 

211 BETA AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Uncertainty in the wind-relative aerodynamic sideslip angle 

212 BETAU ‘^(^WA^ AERODYN Given in the description of AERODYN 

in Section 4.1.2 


Euler angles, yaw pitch (Q-j-), and roll expressed in TOPodetic 

coordinates 

213 TOPEUL(l) OP IP Let [TOP + BOD] = [m.j] which was 

calculated previously for words (101-103) 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


214 TOPEUL(2) 9^ 


215 T0PEUL(3) 


OPIP (1) Pitch not vertical 
i|) = arctan(mj^2/"’ii) 
9 = - arctan(m^2) 


<J) = arctan(m 23 /m 23 ) 


OPIP (2) Pitch vertical 
= 0 

9 = + 90° = - sign(90, m^^) 
(j) = arctan(m2j/m2j^) 

OPIP (3) For Ascent and if ^ < 0°, 

(() = (j) + 360° 


Euler angle rates relative to the TOPodetic frame 


216 

TOPRATE(l) . 


OPIP 

[ECI T0P3 was computed previously for 
words (101-103) 

217 

(2) 

w[T0P] 

OPIP 

[M50 > ECI] was read from input file OPIP. IN 




M3X3 

[M50 > TOP] = [ECI -*• TOP] [M50 ^ ECI] 

218 

(3) ^ 


EULRATE 

u[T0P] = 9^, $y]^ 





calculated by algorithm given in the 


description of EULRATE in Section 4.1.2 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Inertial body angular rates around the BODy axes 


219 

IBRATES(l) X 

OPIP 

Read M50 state vector from input file 
TRJATTDATA: $ = [Si... Sia]^ and 

220 

1 1 

o 

o 

CD 
1 — 1 

(Nl 


convert from radians to degrees 

221 

(3) * 




! “xB ^13 * ^D/R 

V = ^14 * S/R 

“zB " ^15 * S/R 

OP IP m[BOD] = 


Uncertainties inertial body angular rates around the BODy axes 

222 IBRATEU(l) . OPIP Read ^(^[MSO]) = from the input 

I file TRJATTDATA and convert from radians to 

I degrees 

223 (2) > a(i!:[B0D]) 

224 (3) ' 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

IBRATEU continued 

/ ‘^Kb^ " ^^13,13^ ^ * Vr 

OPIP s a(o)^g) = (Cj 4 J 4 ) ^ 

= (^15,15^ ^ * ^D/R 

OPIP o(w[B0D]) = [o(w^g), o(Wyg), 


Inertial body angular accelerations around the BODy axes 


225 

IBANGA(l) V 

OPIP 

Read M50 state vector from input file 


); 


TRJATTDATA; $ = [S^... Sjg]^ 

226 

(2) > o)[B0D] 


convert from radians to degrees 

227 

(3) ' 




! “xB " ^16 * ^D/R 

“yB " ^17 * ^D/R 

“zB " ^18 * S/R 


OPIP 


iS = [(1) 


xB’ '"yB’ *^zB 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 

word name symbol program algorithm 

Uncertainties in the inertial body accelerations around the BODy axes 


228 IBANGAU(l) 


229 


(2) > a(5) 


OPIP Read lt(5[M50]) = [C,j]ig^j 3 

from the input file TRJATTDATA and convert 
from radians to degrees 


230 


(3) 


OPIP 


°(“vr) = (Cl« ^c) 



yB 


D/R 

'17,17^ "D/R 


\l'i 

“zB' " '" 18 , 18 ' * S/R 


(C-ip ip) 


OPIP 


r* • • tT 

P(») = 


Sensed acceleration of the vehicle Center of Mass expressed in BODy axes coordinates 


231 DDBODC(l) 


M3X1 Ig [BOD] = [M50 > BOD] r^ [M50] 


NB 


NB 


232 


(2) > r. [BOD] VADD 


CM 


r.pEBOD] = a [BOD] + a [BOD] 
JJJ CEN ROT 


233 


(3) 


CROSS 


;^[B0D] X Tf^g[B0D] 

CM 

+ J[B0D] X Tf^gEBOD] 
CM 


VADD 


h [BOD] = L [BOD] + TmJBOD] 


CM 


NB 


NB' 

CM 
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BET OUTPUT PRODUCTS (continued) 


DEFINITION OR DESCRIPTION OF ANCILLARY DATA PARAMETER 
word name symbol program algorithm 


Uncertainty in the sensed acceleration of the vehicle Center of Mass in BODy 
coordinates 


234 DDBODCU(l) 


235 

236 


(2) > Mh CBOD]) 


CM 


(3) 


OP IP Read §[M50], r^igCBOO], [M50 > BOD] 

CM 

and $(![M50]) from TRJATTDATA 


OPIP Read a(r|^g[B0D]) from OPIP.IN 
CM 


DDBCU f(r„ [BOD]) = fctn(^[M50], TmoCBOD], 


CM 


NB‘ 

CM 


[M50 ^ BOD], (f(![M50]), a(r..[B0D]) ) 


NB‘ 

CM 


OPIP 


o(rg [BOD]) = C22'^^ 


1 / 2 ^ 


CM 


where [C.^] = [BOD]) 

CM 


Position of the vehicle Center of Mass relative to the Nav Base in BODy coordinates 


237 BODCG(l) 

238 (2) 

(3) 


OPIP 


?nbCbod] 

CM 


Read from the input file TRJATTDATA and 
placed directly into the output buffer 
array. 


239 
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4.1.2 Subroutines Called by OPIP 


This section contains the description of subroutines called by OPIP and which 
contain significant algorithms needed for computing of Ascent/Descent output 
products. The subsection for each subroutine starts with a list of terms in the 
subroutine calling argument. Each output term is indicated in the text by the 
symbol *>> next to the left-hand margin. 

The list of subroutines documented in this section is as follows: 

SUBROUTINE PAGE 


AERODYN 

65 

BANK 

76 

EULANG 

79 

EULRATE 

94 

FPANG 

96 

FPANG2 

99 

GEODGEO 

103 

GEOGEOD 

105 

SPBCALC 

108 

WINDTOP 

no 
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AERODYN 


Calling Argument List (in list order) 


Input 

Li First-pass flag. 


^D/R 


Conversion constant equal to degrees per radian. 


V [BOD] Velocity of vehicle Center of Mass with respect to the 

wind. This vector is expressed in vehicle BODy 
coordinates. 


(f(v [BOD]) 
wi nd 
CM 

ag[B0D] 

CM 

^AM 


“(Tar> 


Covariance matrix of v [BOD]. 

wind 

CM 

Sensed acceleration of the Center of Mass expressed in BODy 
coordinates. 

Atmospheric mass density (slugs/ft^). 

Uncertainty in atmospheric mass density. 

Atmospheric temperature in degrees Rankine. 

Uncertainty in the atmospheric temperature, 


Output 


“W 

^WG 


Wind-relative angle of attack (degrees). 
Uncertainty in angle of attack. 

Wind-relative guidance sideslip angle (degrees). 
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^WA 

q 

o(q) 

M 

a(M) 

''ea 

CM 

CM 


WO 


L 


D 


V 


a(v) 


Uncertainty in 3yjg. 

Wind-relative aerodynamic sideslip angle. 

Uncertainty in 3^^^. 

Total dynamic pressure. 

Uncertainty in the total dynamic pressure. 

Pitch dynamic pressure. 

Yaw dynamic pressure. 

Wind-relative Mach number. 

Uncertainy in M. 

Equivalent air speed of the Center of Mass (knots). 
Uncertainty in 

CM 

Wind-relative load factor (in units of g, the sea-level 
acceleration of gravity). 

Wind-relative drag acceleration (magnitude). 

Wind-relative lift over drag. 

Wind-relative viscous parameter. 

Uncertainty in the wind-relative viscous parameter. 
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The Algorithms 


Constants (set on 
g = 32.174 
= 12.1527 
= 4289.05 
y = [0,1.0]^ 


first pass through the subroutine) 

Earth-surface gravitational acceleration. 

Constant for computing the equivalent air speed, 

CM 

Constant for computing the Mach number, M. 

The (column) unit vector used to compute lift over drag. 


L 


D* 


Preliminary computations: Magnitude of wind-relative velocity and its 

variance 


Let the magnitude of v [BOD] be V; 

wind 

CM 

2 2 2 1/2 
V = (Vx + Vy + V^ ) ' . 

r\W J\W p\U T 

To compute the variance of V, first C-^3 = 3^3 is computed. 

3^ 3Vx avy av^ 


ft 2 2 2 1/2 V 

av \ ^ ''z ^ = IT* 


av 


Likewise, 


av _ ''y av _ ''z 

3v “ V 3v^ ■ V • 


The variance of V is computed from this and the covariance matrix: 


/w \2 _ p 3 V 3 V av -| ||./ rnfirni pay av av -.T 
°vV) ~ Tv”"* [“ODJ) Lay > ciw » ciw J 

X y z 


wind 
CM 


-3V ’ 3v ’ av 
X y z 


The uncertainty, which will be used later, used here is the standard 


67 



deviation given by 


2 1/2 
o{V = Ca(V) ] ' . 


Wind-relative angle of attack, a 


W 


Let 


2 2 1/2 

\z = (% * ) • 


Then 


V V 

, Z X 

sin a, = -r, — and cos a = -n — 
” ''xz “ ''xz 


So, 

*» 


\ = arcsin(-j^) 


xz 


and is converted to degrees and adjusted so that 0 ^ 180°. 


Uncertainty in 


First, the partial derivative matrix is computed: 

r *“« n r’“w ’“U ’“W,T r' “V „ “W-.T 

S; [BOD]-* ■ ^'z-* ■ "iB ’ ’ • 

wind 
CM 


Using this and the covariance matrix, the variance is calculated from 

- r’““ If(„ r«nnn r*°“ 

^ ^av ’ 3v ’ 9V ^ Pv ’ 3v ’ 3V ^ • 

X y z wind x y z 

CM 


The uncertainty, or standard deviation, in is then 


*» »(ok) = 
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Wind-relative aerodynamic side slip angle, 


As before, 

•». 2 
V = |v [BOD] I = (v + V 
wind 
CM 


Then, ^ 

sin and cos = 

So , y 

*» = arcsin(-^) 

and Byjp^ is converted to degrees and adjusted so that 0 ? By^p^< 180°. 
Uncertainty in B^^^^ 


sin sin 

V -■ 


First, the partial derivative matrix is computed: 
r ”*“«A , . r*“«A ’“uA ’"wA,T 

L~i -j - j 


3v [BOD] 
wind 
CM 


■Sv, ■ 3Vj, 


“w ®UA “UA 


2 1 /2 

+ V, ) / . 


2 2 1/2 

(V, . ) / 

^ 


Using this and the covariance matrix, the variance is calculated from 

2 ^^Wa ^^WA, rpQp-,v r^^WA ^‘^WA ^WA ,T 

“^^Wa) = tl77’ winT^^ "''y ’ ^''z 

CM 

The uncertainty, or standard deviation, in B^^^ is then 




WA 


WA^ 
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Wind-relative guidance side slip angle, H 


WG 


2 2 1/2 

Let V =(v +v ) . 

xy ^ X y ^ 


Then , 


V V 

^WG “ ^WG ~ V 

xy xy 


So, 

*» 


^WG " arcsin(/-) 


xy 


and is converted to degrees and adjusted so that 0 < ii^g < 180^ 


Uncertainty in 3 


WG 


First compute the partial derivative matrix 


30, 


80,,- 83,,- 83, T - sin 0,,- cos 3,,„ , 

r WG T _ r WG WG WG-iT _ r WG WG ni' 

L J L^w > 2./ > J - L \# > \# > UJ « 


V- ronm ’ 3v ’ 8v 

8v [BOD] X y Z 

wind 

CM 


V ’ V 
xy xy 


Using this and the covariance matrix, the variance is calculated from 


2 83 ,,- 30 ,,- 33 ,,- ^ 33 ,,- 83 ,,- 33 ,,- t 

fo \ r WG WG WGt rDnnT\ r WG WG WGtT 

“(“wg' ” [ 3 ^-- 3V— > — ] C— ■ — • • 

X y z wind x y z 

CM 


The uncertainty, or standard deviation, in 0^^^ is then 






Dynamic pressures and the uncertainty in q 


The (total) dynamic pressure, q, is calculated by 


*» 


q 2 ^ am '' 
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where is the atmospheric mass density. 

The uncertainty in q is calculated from 

*» »(q) = [p^M^ V »(v)^ ♦ 7 

where is the uncertainty in p^j^. Note that both and '^(p;^h,|} are 

brought in through the argument list. 

The pitch dynamic pressure, q^, is calculated from 

*» q„ = V- 

The yaw dynamic pressure, is calculated frqm 


The equivalent air speed, V-., and its uncertainty, <5(Vr«) 


'EA 

CM 


EA^ 

CM 


The equivalent air speed, in knots, of the Center of Mass is calculated by 

''ea " *^ew ‘^AM ^ 

CM 




V. 


Its uncertainty is calculated from 
*» 




CM 


The load factor, n 

The load factor is computed by dividing the sensed acceleration of the Center 
of Mass, a^CBOD] brought in through the calling argument, by the Earth's 
CM 

(surface) gravitational acceleration, g. Thus, 
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*» 


lag)[BOm| 

CM 


The wind-relative drag acceleration, 

First, the unit vector of the wind-relative velocity of the Center of Mass is 

calculated; let it be v . Then 

w 

V [BOD] 


wind 

CM 



And then the magnitude of the drag acceleration is computed from 
^WD 

CM 

The wind-relative lift-over-drag, Lp 

»-- - ■■ - ■ — — , . . ■ — ■ - . . . . .1 , — i— 

^ T 

Lift is computed as follows, y = [0,1,0] was set up during the first pass 
through the program along with other constants. The direction, L, of the lift 
force is determined from the cross product 

->■ 

L = y X v^j. 


The unit vector in this direction Is 



The lift acceleration-, a^, is the component of a^gCBOn] along L; thus, 

CM 

a^ = a^EBOn] • L. 

CM 

The lift-over-drag is the ratio of a|_ to if a^jp * 0. Thus, 

^ ^WD ^ 

“ " "hd ' ° 
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Mach number, M, viscous parameter, v, and their uncertainties 


First, the ambient atmospheric temperature and its uncertainty is converted 
from degrees Rankine to degrees Kelvin; 

r = — T = T 
'ak 9 AR a 

and 

Next, the stagnation temperature, T^, and its uncertainty are computed using 
V, the magnitude of the wind-relative Center of Mass velocity. The algorithm 
i s 

b 2 

= 726.97 + 0.468 + 3.4098447*10" V . 

Its lincertainty is computed from 

a(T^) = [(0.468 a(T^) + (4.650816351*10"^^ V o(V) 


Next, the "intermediate coefficient of viscosity", C^, and its 
uncertainty, o(C^), are calculated using the following algorithms: 


C 


0 


Tc 1/2 

(^) ' [T, + 122.1(10 '')] 

'a ^ 

[Tg + 122.1(10"^^^$)] 


and its uncertainty, 


T <j(T >/2 

»(C„) - [-^ + ^ ^ -] • 

'‘^A 


Next, the coefficient of dynamic viscosity, u, and its uncertainty, o(m), are 
calculated using the following algorithms based on the U. S. Standard 
Atmosphere of 1976; 
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T, 


3/2 


„ • 3.0449939‘in-“ (rpr^jj^) 


and its uncertainty, 

8 1 /2 {-j + 165.6) 

o(m) = 3.0449939*10 T. ' a(Tj — ^ r- 

^ + 110.4) 


Finally, the wind-relative Mach number and its uncertainty are computed using 
the following algorithms: 




M = V(k, T,) 






and its uncertainty, 
a(v) 


o(M) = [ 


2 a(T,)" 1/2 

+ / 


17156.2 t/ 


where, as before, V is the magnitude of Vg[B0D] and o(V) is its uncertainty. 

CM 

The constant kj^ is set during the first pass through the routine. 


The last output terms to be computed are the wind-relative "viscous 
parameter", v, and its uncertainty, o{\>). They are computed with the 
following algorithms based on the U. S. Standard Atmosphere of 1976: 




V = M[- 


^0 ^ 


1 /2 


107.5 V p 


AM 


Its uncertainty is computed from the following series. 


These terms are computed; 


a(M) Cq p 


'pss 


107.5 V p 


AM 
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'pss 


(107.5 Vp^^) (4 Cq) 


2 2 
a(p) M C 


0 


(107.5 V P^^) (4 p) 


V 


2 2 

V M Cq p 


PSS 3 

(107.5 (4 V ) 


P 


pss 3 

(107.5 V) (4 p ) 


With these, the uncertainty in the "viscous parameter" is 


1 ,2 

*» a(v) = (M + C + P + V + P ) ' , 

^ ^ pss pss pss pss pss' 
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BANK 


Calling Argument List (in list order) 
Input 


^'n/R 

'-ASC 

[TOP > BOn] 

V [TOP] 
wind 
CM 


Conversion constant equal to degrees per radian. 

Flag which is .TRUE, if ascent analysis is being done. 

Transformation matrix, TOPodetic to BODy frame. 

Wind-relative velocity of the Center of Mass expressed i 
the topodetic frame. 


Output 

$ Wind-relative bank angle in degrees. 

w 

The Algorithm 
First, the unit vector 

V [TOP] 
wind 

" = CM 

|v [TOP] I 
wind 
CM 


is calculated. Next, the body pitch and roll axes are extracted from the 
transformation matrix. Let 
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[TOP ^ BOD] = 



This means that each row of the TOP to BOD matrix is the row unit vector 
representation of a rotation axis expressed in the topodetic frame. Expressed 
as the column vectors, they are 

r = roll body axis expressed in the topodetic frame; 
p = pitch body axis expressed in the topodetic frame; and 
y = yaw body axis expressed in the topodetic frame. 

y\, /V ^ 

Let n be the unit vector normal to both r and v^. The normal vector, n is 
found from the cross product, 

n = r X V 

Its unit vector, which has the components nj, n2, n3, is 

^ n T 

n = — = [ni, n2, ns] . 

|n| 

^ A 

The angle, a, between r and v^^ is found from 

A* 

a = arccos(r . v^). 

A ^ 

Then the rotation matrix, [R] which rotates r about the axis n through the 
angle a is: 


[R] 


[ni + (1-ni) cos a] [nin2(l-cos a) - 03 sin a] [nin3(l-cos a) + n2 sin u 

2 2 

[nin2(l-cos a) + n3 sin a] [n2 + (l-n2) cos a] [n2n3(l-cos a) - n^ sin a 

2 2 

[nin3(l-cos a) - n2 sin a] [n2n3(l-cos a) + ni sin a] [n3 + (l-n3) cos u] 
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Next, determine the new pitch axis, P, after the rotation. It is found by the 
following logic; 


If V [TOP] does NOT lie along r; P = [R] p 
wind 
CM 


If V [TOP] DOES lie along r; P = p. 
wind 
CM 

Define w as the unit vector in the topodetic x-y plane which is also normal 
to Vy. If is expressed in terms of its components, 

then 

w = [-Vy, v^, 0]^ 


and its unit vector is 


Finally, the wind-relative bank angle, is calculated by 
*» = |arccos(w • P)|. 


This, being in radians, must be converted to degrees. 
Finally, adjustments must be made. Given that P = [P^, P^, 


if P^ < 0, change the sign of 

Also, if this is an ascent analysis (L^ 5 q = .TRUE.) and < 0, then add 360° 
to $ 1 ^. 
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EULANG 


Calling Argument List (in list order) 


Input 


S/R 

'-ASC 




Conversion constant equal to degrees per radian. 

Flag which is .TRUE, if ascent analysis is being done. 


Flag which is .TRUE, for the first pass through the routine. 

Error counter (This can be incremented within the routine, 
so it also can become an output). 


The 15 -element state vector 

[Si. S2, Saf = fg[B 0 D] 

NB 


[S4, S5, Ssf = V [BOD] 
wind 
NB 


[$i ••• Sis]^. 

Geocentric position vector of 
the Nav Base in body 
coordinates. 

Wind-relative velocity of the 
Nav Base in body coordinates. 


[S?* Se, 



= a [BOD] 
wind 
NB 


Wind-relative acceleration of 
the Nav Base in body 
coordinates. 


Sio» Sii, S12 


Not used. 


${S) 


[$13, $14, $15]^ = w [ECI] Inertial (ECI) angular rates 

D 

of the body axes. 

Covariance matrix of the state vector, §. 
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w[BOD] Wind velocity expressed in body axis coordinates. 

(j;(w[B0D]) Covariance matrix of v [BOD]. 

W 


Output 

'^W 

<‘u 

♦w 

♦n 

“(V 

<=(♦„) 


Wind-relative yaw angle in degrees. 
Wind-relative pitch angle in degrees. 
Wind-relative roll angle in degrees. 
Uncertainty (standard deviation) in 

Uncertainty (standard deviation) in 9^^. 

Uncertainty (standard deviation) in (jiyj. 

Wind-relative yaw rate in degrees/sec. 
Wind-relative pitch rate in degrees/sec. 
Wind-relative roll rate in degrees/sec. 
Uncertainty (standard deviation) in 

Uncertainty (standard deviation) in 

Uncertainty (standard deviation) in 
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Algorithms 


Determination of a wind-related coordinate system 


First, obtain the relative velocity v [BOD] = v^ which is 

REL 


; ■ 5 [BOO] - S[B0D] = [V V V f 

Wind ^ 

NB 

^ ^ ^ 

Then, define a coordinate system x, y, z such that 

y^^ is normal to Vj^ and r^CBOD]: the "right" vector; 

NB 

is opposite in direction to r^CBOD]: the "down" vector; 

NB 

and is normal to y^j and z^^: the "forward" vector (which completes 

the triad). 

For notational simplicity, let r^[B 0 D] be represented by fg = [xg, yg, Zg]^. 

NB 

Since uncertainties of the wind-relative Euler angles are to be calculated, it 
is necessary to calculate partial derivative matrices. They are set up in 
subroutines P 3 UV and PCROSS. The following notation will be used. 

If a = [ai, a2, 83]^ and S = [bi , b2, b3]^ are two vectors, then the partial 
derivative matrices will be indicated by 


3ai 

3ai 

3d 1 

3 b ^ 


W 

3a2 

3a2 

3d2 

3b 

357 

W7 

3aa 

3aa 

3as 

3F7 

357 

3F7 
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Compute the "right" unit vector, y, and its partial derivative matrices 


The subroutine PCROSS is called to compute the cross product 


y = V. X f„ = V [BOD] X rjBOD]. 


REL 


NB 


Then it computes the skew-symmetric partial derivative matrices of y with 
respect to the two input vectors, Vj^ and r^, of which y is the cross product, 


[i^] and 


3y 


9Vp 


ar 


B 


Another subroutine, P3UV, is called which computes the unit vector of 

y = = Lvi. y 2 . yal 

|y| 

and the partial derivative matrix of y with respect to y: 

9y 

Then, the following partial derivative matrices are computed: 


= [^3 [^-3 


ay ar 


B 


and 


[P(^-)3 = [|^3 C-^3. 

Vn ^y 9Vp 


Compute the “down" unit vector, z, and its partial derivative matrix 


The subroutine P3UV is called to calculate 


+ 

rr 


3r, 




9rr 
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Since the "down" direction is opposite to that of fn, 

D 

^ ^ T 

Z = - Tg = [Zi, Z2, Z3] 

and its partial derivative matrix is 

7 

[P(l-)] = - 

'"b *''b 

A 

Compute the "forward" unit vector, x, and its partial derivative matrices 

Since the "forward" vector x completes the triad, the subroutine PCROSS is 
called to calculate 

AAA 

X = y X z = [xi, X2, X3] , 

[p(4)] = C^] 

and y 

A A 

[p(4)3 = C^3. 

Z 3 Z 


Compute the wind-relative pitch angle, and its partial derivatives 


*» 


®W ■ - '^D/R ( 


2 2 1/2 
(Xl + Zi ) 


). 


0,, is constrained so that -180 ? e,, ^ 180. 
W w 


2 1/2 

(xi + Zi ) ' 


0 , 0 ] 


The partial derivative matrices are expressed as row matrices: 

30 30,, 30,, 30,, Zi 

r_w-| = r_Jl' _W = r L 

Lg- J 8 X 2 ’ 3X3-* •- 

and 

30,, 30,, 30,, 30,, - Xi 

r—Wi = r— ^ = r 

^ J •-3Zi’ 3Z2* 3Z3-' ^ ^ 


2 2 1/2 
(Xl + Zi ) ' 


0 , 0 ]. 
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These will be used to determine the angle rates and uncertainties. 


Compute the wind-relative yaw angle, and its partial derivatives 


'•'W " ^D/R (yi) 


Its partial derivatives with respect to y are 

3y 


r 'w, . = r 

L J La., » aw„J ^ L~ 


1 


■9yi’ 9y2* 9y3‘ 


z 2 1/2 

(y 2 + ys ) 


, 0 , 0 ] 


Compute the wind-relative roll angle, <J>^, and its partial derivatives 


*» 


♦w ' - W I 


ya 


2 2 1/2 

[yz + ya ) 


0 


is constrained so that -180° ^ ^ 180°. 


For ascent and <j>^^ < 0, is made positive: = 'i’y + 360°. 


Its partial derivatives are 

3 (|) 

9y 


r W-. 

L - J - l-ayi’ 3y2’ Syg-* 


= [ 0 . 


-yz 


ya 


2 2 2 2 

yz + ya yz + ya 




Compute partial derivative matrices used to compute angle rates and their 
uncertai nties 


The matrix of the partial derivatives of the wind-relative Euler 

A ^ 

angles, , with respect to the unit vectors y and z is calculated to be 
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3E,,. 

3(y, z) 


% 


3»U 

3'1'w 



■5yT 

3y2 

9Y3 

3Zi 

3 Z 2 

3Z3 

36.. 

90.. 

3 0., 

96.. 

3 0., 

9 0.. 

W 

w 

W 

W 

W 

W 

3yi 

3y2 

9y3 

3Zi 

3^2 

3Z3 






“♦w 

3yi 

9Y2 

9y3 

3Zi 

3Z2 

3Z3 


The individual derivatives, or elements, of this matrix are as follows: 


9'1'w _ 


1 


3»H _ 


- 

9yi 

2 

(yi 

2 1 / 2 * 
+ Y 2 ) 

9y2 9y3 

3Zi " 

3Z2 

9Z3 

9Yi = 

3xi 

3Xi 36^j 

3y . 3X2 

3X2 90y 

9y.j ^ 9 x 3 

3X3 

9y/ 

i = 

1,2,3 

'®W 

^^i 

3Xi 

3Xi 36,1 

+ ^ 
3z^ 3X2 

3X2 90^ 

3z. 3x7 

3X3 

i = 

1,2,3 

_ 

0; 

94'w - ^2 

!!w = _ 

Y3 



9yi 

9yz 2 

2 ’ 9 y 3 

2 

2 



= 0 


Y2 + Y3 


Y2 + Y3 


and 




__w 

3Zi 


0 

aZ2 3Z3 


The matrices of partial derivatives of the state vector position, and 

wind-relative velocity, Vn, with respect to the wind-relative Euler angles, 

^ T 

E^j^, are as follows. With rg = [Xg, yg, Zg] , 
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Another matrix that is used for the determination of Euler angle rates is that 

containing the partial derivatives of y and z with respect 

to Og. Vg, ♦„. WUh fg . [Xg, yg, Zgf and v = [Vg^, Vg^. »gpT_ 


-9(y» z) 


3 (^ 8 * Vj^, iPy, 0^, 


3 = [P,J 
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9yi 

9yz 

3y3 



3Xb 



9yi 

3y2 

9ya 




’>'b 


3yi 

ay2 

ays 

[p,j] - 


'^B 

'^B 


3Zi 

3Z2 

3Z3 



’"b 

33^ 


3Zi 

3Z2 

3Z3 



'^B 

^^B 


3Zi 

3Z2 

3Z3 



'^B 


The elements 

Pij ’ 

= 1 

- 6 . 

Let 

’ ^2 = 


and 

Pi. 

3 

. = M 
J k=l 


3r 


Pl7 Pl8 Pl3 


P27 P26 P23 


Pa/ P38 P39 


0 0 Pit/ P48 P 49 


0 Pb/ Pb8 P 59 


0 0 0 P67 P68 P69 


-^3 


Wk ^ 

^ w ' 


and for i = 4,5,6 with j = 7,8,9 
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3 3z. 

Pii = ^ (tt 


3E 


Wk 




). 


= V 


The index k is used as follows: for k=l, P|^ = Xg and v^, 

Tk = yg and = Vgy; and for k=3, r^ = Zg and 

The next step is to calculate the partial derivative matrix 


Rx 


C-7 


3(rg, V|^, ipy, (j.^) 


-] = ^W> 


3(y. z) 


’R'' 


= tPij] 


where 




Pi 2 


Pi3 


Pi^ 


Pi 5 


Pi6 



3yi 


*Ewi 

dZi 



+ 

A 


9y 



3Z 


'^Wi 

syz 


^^Wi 

3Z2 

— X 


+ 

A 


3y 



3Z 

^^B 

'^Wi 

3ya 



3Z3 

A 


+ 

A 

T“» 

3y 

^"b 


3Z 

'^B 

^Ewi 

3yi 

+ 

’^Wi 

3Xi 


X" 

3y 


+ — X 

3Z 

- A 

ayz 


3y 


3Z 

— A 

3ya 


3y 


3Z 

elements 

are 

taken ' 


3Z2 

“T“» 


3z 


B 


3v„ 


1J- 


Pij " Pij ^ ^ ^ "" 7,8,9. 


; for k=2. 


3(y». i). ] 

’W'*' ’W^ ’W*** ^ 
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The angle rates, 0^^, and are calculated from 


u»[B0D] = 


®W’ 


T 


= [• 


®W’ 


3(rg, v^, ♦„) 


] “z’ 'x- 




Note the right hand matrix is a portion of the state vector, §, which includes 
words 4 through 9 and 13 through 15. 


Specifically in terms of the elements Pij 


CPii3 = C— . 

9(rg, V|^, «p^, 4>y) 


]: 


the angle rates are 

*» = PiiVj^ + Pi2Vy + Piav^ + + Pisv^ + Piev^ + Pi7tj + Pie^j + Pisi*! 


*» = P 2 lVj^ + P 22 Vy + P 23 V 2 + P24Vj^ + P 2 SVy + P 26 V 2 + P27'i'j + P2b‘^i + P2bVj 

*» <}>w " ■*■ ’'32Vy + P 33 V 2 + P34Vj^ + PasVy + P 36 V 2 + P37'l'i + P38®i + P33<i>i 

These are converted from radians per sec to degrees per sec. 


Compute the wind-relative Euler angle rate uncertainties 

Consider the [P] matrix as a three-element (column) matrix whose elements are 
nine-element row vectors; that is. 


89 


[- 


*w> ♦w) 


'r’ *W’ °w 


-] = [i*i, h, fsf 


where 


P. 


•j “ CP-ji> 


The state vector is dimensioned as a 15-element array, and, therefore, its 
covariance matrix, (|)(S), is dimensioned 15 x 15. However, elements 10, 11, 
and 12 are not used, so, for the Euler angle and angle rate uncertainties, 
terms of $ and $($) involving position, Tg, are not used. The elements 
of ^ and (p(?) are accessed in such a way as to effectively make § a nine 
element vector and (f($) a 9 x 9 matrix. 


Then consider (|:($) as a nine element column matrix of nine element row 
vectors, 

(K^) = [Cij] = [Cl, C2 ...Cgf = [5j], 

where 

~ [C j 1 , ...Cjg3. 

Now, define as the dot product 

R. .=?-•* 
ij J 1 

Nine of these, 1 ^ j ^ 9, are arranged to form a vector, . 

The variance, of the i^*^ wind-relative Euler angle rate is computed 

from 


= Pi • Pi-, i = 1.2,3. 


In other words. 


o('l'lj) = Pi • Ri; o( 9 u) = ^2 • R2; a( 4 >y) = P3 * R3. 
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All three variances are tested for being positive. Any which is found to be 
negative is set to zero, and a message is written to an output file. 


Finally, the uncertainties (standard deviations) of the wind-relative Euler 
angle rates are computed from 

*» ^ ; 

*» 0(0^) = and 

*» a(4>w) = 

These are all converted from radian/sec to degrees sec. 

Compute the wind-relative Euler angle uncertainties 

Again, considering (j)($) as a 9 X 9 matrix, the wind velocity covariance 

matrix, (f(v [BOD]), which is a 3 X 3 matrix, is added to it. If 
wind 

${v [BOD]) = [v^..], 
wi nd 

the sum, (|)(i!i), a 9 X 9 matrix, is formed as follows: 

It(fi) • »(§) ♦ ((v [BOO]) 
wind 
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m = 


C 34 


C35 


C36 



+ 

Vll 

Ci *5 + 

V12 

C 46 + 

Vi 3 

• • • 

Cs 4 + 

V 21 

C55 + 

V22 

C56 + 

V23 

• • • 

Cet* + 

V3I 

Ces + 

V32 

Cee + 

V33 

• • • 

C?*t 


C75 


C?6 



• 


• 


• 



• 

• 


• 

• 


• 

• 


C95 


(I)(W) is considered as a column vector with each of the nine elements being a 
nine-element vector. Let 


If(ij) = Ci5i. 
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All three variances are tested for being positive. Any which is found to be 
negative is set to zero, and a message is written to the output file. 

Finally, the uncertainties (standard deviations) of the wind-relative Euler 
angles are computed from 

*» 

*» cr(0w) = and 

*» a(4.y) = 

These are all converted from radians to degrees. 
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EULRATE 


Caling Argument List (in list order) 
Input 


rg[M50] 

NB 

rg[M50] 

NB 

iUCBOD] 

[M50 B0D3 

[M50 T0P3 
[TOP B0D3 


Geocentric position of the Navigation Base in M50 
coordinates. 

Geocentric velocity of the Navigation Base in M50 
coordinates. 

Orbiter angular velocity in body axis coordinates. 

Transformation matrix from M50 to BODy axis coordinates 
(The Navigation Base is the origin of the BODy system.). 

Transformation matrix from M50 to TOPodetic coordinates. 

Transformation matrix from TOPodetic to BODy coordinates. 


Output 

w[T0P] 


Topodetic Euler angle rates for yaw, pitch, and roll: 
u[T0P] = [i|>y, 6^, 


Input 

Conversion constant equal to degrees per radian. 
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Algorithm 


The algorithm is described in TRW IOC 83:W482.4-28, dated 25 March 1983 and 
written by Darwin H. Poritz. This document is included in the Attachment. 
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FPANG 


Calling Argument List (in list order) 
Input 


^D/R 

V [TOP] 
wind 
CM 


o(v [TOP]) 
wind 
CM 


Output 


a(Y„) 


Conversion constant equal to degrees per radian. 

Wind-relative velocity of the Center of Mass in the 
TOPodetic frame. 

Uncertainty in v [TOP]. 

wind 

CM 


Magnitude of v [TOP]. 

wind 

CM 

Uncertainty in 

Wind-relative ("local") flight path angle in degrees. 
Uncertainty in 

Wind-relative ("local") azimuth, or heading angle, in 
degrees. 




Uncertainty in 
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A1 gorithms 


Compute the magnitude, and its uncertainty, 

CM 


Then, the magnitude of this velocity is 


*» 


2 2 2 1/2 
''WR " ^''Rx ''Ry ''d^ ) 


'Rz 



Then, the uncertainty, or standard deviation, in V^j^, the magnitude of the 
wind-relative velocity, is 


*» 


o(Vkr) = [■ 


“Rx “v(v 


x'WR 


2 

) + V 


Ry_^' 

'^WR 


UR 


2 

) + V 


Rz 




WR^_jV^ 


Compute the wind-relative flight path angle and its uncertainty 

The wind-relative flight path angle is computed from 

*» Yy = arcsin(— ^). 

“ ’WR 

which is converted from radians to degrees. 

Its uncertainty, or standard deviation, is computed from 




*» 


222 222 22 

• *^x^^WR^ ^Rx ^Ry ^ °y^%^ ^Ry ^Rz °z^^WR^ ^^Rx ^Ry ^ 


'^WR ^''Rx ''Ry ^ 
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Compute the aximuth, or heading, and its uncertainty 


The wind-relative aximuth angle is computed from 

Vr 


*» 


= arctan(-^) 
''rx 


which is converted from radians to degrees. Its uncertainty, or standard 
deviation, is computed with 

*» oi'V ) = ^Ry ^ °y(^WR) ^Rx ^ ' 


2 2 

''r 

X y 
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FPANG2 


Calling Argument List (in list order) 


Input 

^D/R 

r®[M50] 

NB 


Conversion constant equal to degrees per radian. 

Geocentric position of the Navigation Base in the M50 
frame. 


r^[M50] Geocentric velocity of the Navigation Base in the M50 

NB frame. 


Output 

R|^g Magnitude of fg[M50]. 

NB 

Magnitude of r^CMSO]. 

NB 

Yj Inertial flight path angle measured, in degrees, in 

horizon plane 


Inertial azimuth, or heading angle, measured in degrees, 
from the North. 


Algorithms 

Computation of the magnitudes of position and velocity 

Let fgCMBO] = [Xj, yj. z^f and rg[M50] = [v^, v^. v^]^. 
NB NB 
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Then, 


2 1 


*» 

*^NB ~ 

(Xj + Yj + ) 

and 





2 2 2 1 

*» 

'^NB " 

+ “y * ^ ) 

Computation of 

the inertial flight 


/2 


First, the unit vectors r^^g and Vj^g for the position and velocity vectors, 
respectively, are computed: 

^ ^IT 

''nB ^R|^g’ Rj^g* Rjjg^ 


and 


'NB 


- -^1 • 
''nb ''nb ''nb 


The inertial flight path angle is computed from 

*» (Sb • 'nb> 

which is converted from radians to degrees. 


Computation of the inertial azimuth or heading, 


First, the cross product of velocity and position is calculated; 

t = rg[M50] X rgCMBO]. 

NB NB 


Then, the horizontal velocity, Vg, is calculated from 

Vg = r^CMBO] X t 
NB 
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and its unit vector, v^, is determined; 



The vector, K, pointing North is computed from 

= rg[M50] X (k X rg[M50]) 

NB NB 

where k is the (unit) vector [0,0,1]^, 

Then, 

A 

N = (Xji + Yjj + Zjk) X [k X (Xji + y^j + Zjk)] 

A ^ 2 2 

= - zx i - zyj + (x + y )k 

2 2 T 

= [- zx, - zy, (x + y )] . 

The unit vector is calculated from 



Next, the unit vector, E* pointing East is computed from 
E = il X 

Then, 

A A /V A 

COS '•'i " N • v^ and sin Tj = E • v^. 

To avoid a singularity when computing fj, cos fj is checked. 
If cos 'Pj = 0 and sin 'Fj > 0, set fj = 0° 
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If cos = 0 and sin ?0, set "Jj = 270° 
If cos 0, then f is computed from 

sin f, 

*» ’l • 

which is converted from radians to degrees. 
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GEODGEO 


Calling Argument List (in list order) 
Input 



Conversion constant equal to degrees per radian. 

Array of Earth geophysical parameters. Those used here are 
as follows: 

Rgg = Earth's equitorial radius = Pg{l) 

Rgp = Earth's polar radius = P£(6) 

Geodetic latitude, in degrees, of the desired location, P. 

Geodetic longitude, in degrees, of the desired location, P. 

Geodetic height, in feet, of the desired location, P. 


Output 

fg[GE0] Geocentric position vector of the desired location, P, in 

P the GEOgraphic ("Earth-Fixed-Greenwich") coordinate system. 


A1 gorithm 

First «|)p and are converted from degrees to radians. Then 


and 


j _ - ^9E 

xy RIJ TTii 

[1 + tan 4 .^) ] ' 
®E 
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^®P 2 

*^xy *D’ 


Then, the geocentric position vector is 


rg[GE0] 

P 



(dxy + hj COS ♦(,) cos 
<“xy * % *0> 
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GEOGEOD 


Calling Argument List (in list order) 


Input 


C 


D/R 


R 


OE 


'^«P 

rgCGEO] 

P 


Conversion constant equal to degrees per radian. 
Equitorial radius of the Earth. 

Polar radius of the Earth. 

Position vector of a point, P, (eg., for the vehicle, the 
Center of Mass, or the Navigation Base) expressed in 
GEOgraphic coordinates. 


Outputs 


X 

h 

5 


D 

D 

D 


VflWV/VIC i* I V. 


1 

I U U I VUVIC ) 


!•>% ^ 

I II ViC^t CCO 9 


rs ^ 4 

Ul pu I IIU 


n 

r • 


Geodetic longitude, in degrees, of point P. 
Geodetic height, in feet, of point P. 
Declination, in degrees, of point P. 


Algorithms 

Preliminary calculations 

First, the Earth's flattening, f, is calculated from its definition: 
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The components of the geographic position vector are Xg, y^, z^; that is, 

TgCGEO] = [Xg, yg, Zg] . 

P 


The following parameters are calculated from f and the geocentric position: 

2 2 
A = Xg + yg 

and 

2 2 

D = (1 - f) Zg . 

Now, a parameter, Bp^, is calculated by means of a five-pass iteration. The 
iteration is initialized by setting 

B = 0.0067 
0 

Then, for i = 0,1, 2, 3, and 4 


( 2 ) 

(3) 


- f) R,e 

i 1 /2 

(C + D) / 


Steps 1 through 3 are repeated until i = 4. 
Then, 
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Calculate the geodetic latitude, <j>p, and longitude, 


The latitude is calculated from 


*» <|)q = arctan[- 


2 1/2 
(A/Bpi) ' 


] 


which is converted from radans to degrees, 


The longitude is 


*» X- = arctan(-^) 

U Xg 

which is converted from radians to degrees. 


Calculalate the geodetic height, h^ 


The height is calculated from 




B.n - fl 

h _ ■ m - f) ^ 


Calculate the declination, 6 


The declination is calculated from 


*» 


6 = a retail [• 


2 2 1/2 

(v, ) 


] 


which is converted from radians to degrees. 
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SPBCALC 


Calling Argument List (in list order) 

Input 

(l:(e[B0D]) Covariance matrix of errors in body axis orientation which 

are errors in the body axis Euler angles in radians. 

e[B0D] = [e^, Eg, where = yaw, e = pitch, and = roll. 

u Vector for which the covariance matrix is to be corrected 

for body axis orientation errors. The vector has the 
elements u = [ui, U 2 , 03 ]^. 


Input and Output 

(P(u) Covariance matrix of the vector u. It is input without the 

effect of body axis errors and output with these effects 
added to it. 


Algorithm 

First, the partial derivative matrix is constructed as follows: 

0 -U3 U 2 

U3 0 -ui 
-U2 Ui 0 

Then, the contribution of the axis uncertainty is calculated by 

3u- 3u- f 

«(=[!]) = [- 57 ^] lf(=[B0D]) [^3^ 
j j 

The axis uncertainty matrix is then added to (f(u) to obtain the output; 


3u. 

J 
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*» ({(u) = 
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WINDTOP 


Calling Argument List 
Input 


^D/R 

Wh 

o{Wn) 

a(6H) 

0(Wy) 


Conversion constant equal to degrees per radian. 
Horizontal wind speed. 

Direction angle of horizontal wind velocity in degrees. 
Uncertainty in W^, 

Uncertainty in 0^ {0^ in degrees). 

Uncertainty in vertical wind speed. 


Hut nut 

~ ^ r 


V [TOP] 
wind 

$(v [TOP]) 
wind 


Mind velocity expressed in TOPodetic coordinates. 

Covariance matrix of v [TOP]. 

wind 


Algorithm 

Calculation of wind velocity expressed in topodetic coordinates 
First, and °(®^) are converted to radians by dividing by Cpyi^. Then 

*» V [TOP] = 

wind 


-Mj^ cos 
-W^ sin 0 ^ . 

0 
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Calculation of the covariance matrix of v [TOP] 

wind 

The matrix elements are set as follows: 

*» [TOP]) = 

wind 

{[o(w^) cos 0^] + [W^ o(e^) sin 0^] } {cos 0^^ sin ° 

(cos 0^ sin 6j^[o(W^) + a(0j^)^]l |[a(W^) sin 0^]^ + [W^ o(b|^) cos 0^]^} 0 

0 0 
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4.2 THE NAVBLK FILE 


The purpose of the NAVBLK file Is to provide the user with a magnetic tape 
file which contains a list of input parameter values used by the trajectory 
estimation (BET) system. This file is produced by the program I2IIPC. The 
contents of the NAVBLK file are gathered by the program INPUT whose inputs are 
mainly user inputs. The only calculated parameters are the declination 
(DELTA) and geocentric distance (RSUBO) of the Nav Base as the vehicle sits on 
the launch pad before launch. These two terms, DELTA and RSUBO, are 
calculated in INPUT; the algorithms for them are given in Section 4.2.2 of 
this document. 


4.2.1 Definition of Terms 

The NAVBLK file contains constants required for a BET solution. The terms in 
this file are defined in Table 1 on the next page. Section 4.2.3 shows how 
the contents of the NAVBLK are arranged in the file records. 
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Term Name Description Source 

IHEADR [180] NAVBLocK header from !»1ET tape See CMET[20] 
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4.2.2 Computation of the Terms DELTA and RSUBO 

These terms are computed for ascent analysis in the program INPUT. RSUBO is 
the geocentric distance of the Navigation Base as the Shuttle sits on the 
launch pad. DELTA is the declination (or geocentric latitude) of the 
Navigation Base before launch. The following symbols will be used: 

Equitorial radius of the Earth. 

Polar radius of the Earth. 

Geodetic latitude of the launch site, in degrees. 

Geodetic height of the Nav Base, at the launch site, above 
the Fischer ellipsoid. 

r RSUBO 

0 

6 DELTA 
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V ^ 

Z 2 


tan <t> 


R,pD + tan ♦ )']*/' 


Let 


/ 2 2 1/2 

V MVxy - ^ 


Then, let 


xy 


and 


V * 


Then, the declination is 


*» ^ " arctan( - / ° ). 

xy xy 0 


The geocentric distance of the Nav Base is 


*» 


.2 21/2 


= C(s, + U ,X„) + (S + u X ) ] 
0 ' z 2 O' ' xy xy o' 
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4.2.3 NAVBLK File Format 


Each record written by the program I2UPC is 1000 words long. The file format 
is shown in the following table. 


Table 2. NAVBLK File Format 


START END 

TERM RECORD WORD RECORD WORD 


IHEADR 1 

GRR 2 

CGT 2 

DPSET 2 

REF 6 

TRAC 6 

CMET 6 

KAPPA 6 

PHIO 6 

DELTA 6 

LAMO 6 

RSUBO 6 

XSUBO 6 

NAVO 6 


1 

1 

5 

805 

5 

32 

632 

632 

633 

634 

635 

636 

637 

638 


1 

2 

2 

6 

6 

6 

26 

6 

6 

6 

6 

6 

6 

6 


180 

4 

804 

4 

31 

631 

631 

632 

633 

634 

635 

636 

637 
640 
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APPENDIX A 

TOPODETIC EULER ANGLE RATES 
(Applicable Document 5, Section 3.0) 


A-1 



83:W482.4-28 
25 March 1983 
Attachment, Page 2 


1. INTRODUCTION 

The Ascent/Descent Best Estimated Trajectory (BET) which is produced by 
the OPIP program provides space in words 216, 217, and 218 of each record for 
rates of Orbiter Euler angles relative to the instantaneous topodetic axes. 
This memorandum proposes a method for computing these angular rates and 
presents this method in Section 3 below. The definition of the Euler angles 
is given in Section 2. 


2. TOPODETIC EULER ANGLES: YAW, PITCH, AND ROLL 


The topodetic coordinate system is defined and displayed in a figure in 
Appendix A of Reference 1. The origin of the system is the point of interest, 
e.g., the Orbiter 's center of gravity or navigation base. The unit vector 
f ^2 representing the topodetic z axis points down (in the direction of 
Earth's gravity) from the origin along the local perpendicular to the Fischer- 
ellipsoid model of the Earth's surface. The unit vector _fj representing the 
topodetic X axis is perpendicular to and points northward in the origin's 
meridian plane. The unit vector representing the topodetic y axis 


completes the nght-handed, orthogonal 


SysLem. 


The Orbiter's body-axis coordinate system is also defined and displayed 
in a figure in Appendix A of Reference 1. The origin of the system is the 
Orbiter's center of gravity. The unit vector representing the body x 

axis is parallel to the Orbiter structural body axis and positive toward 

the Orbiter's nose. The unit vector representing the body z axis is 
parallel to the Orbiter's plane of symmetry, is perpendicular to , and 
points downward with respect to the Orbiter's fuselage. The unit vector 
representing the body y axis completes the right-handed, orthogonal system. 


The topodetic Euler angles of yaw \|i , pitch 0 , and roll <)) are 
defined by fictitiously aligning the Orbiter's body axes x , y , and z 
with the topodetic unit vectors _f ^ , jf^ » ^od * respectively, and then 

rotating the Orbiter in a particular way until the body axes x , y , and 
z are aligned with the body-axis unit vectors » ^rd > 

respectively. The first rotation is about the fictitiously aligned body z 
axis through the angle . The second rotation is through the angle 9 
about the fictitiously aligned by y axis after the first rotation through 
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t|) . The third rotation is through the angle (j) about the fictitiously 
aligned x axis after the rotations through and e . 

Let , J(e) , and K(i()) denote the orthogonal transformations for 

the successive rotations of yaw, pitch, and roll by ip , 9 , and radians 

respectively. From simple figure drawings, it can be shown that 



COS\|» 

si nip 

o" 

I(t) = 

-si nip 

COSip 

0 


0 

0 

1 . 


" cose 

0 

-sine 

J(e) = 

0 

1 

0 


sine 

0 

cose 


" 1 

0 

0 ■ 

K(4>) = 

0 

j 

COSip 

sin^t 


lO - 

si nip 

COSip _ 


The complete rotation L(^,e,())) is given by the product 


L(^,6,(j>) = K{<t>) J(0) I(\|») . 


3. TOPODETIC EULER ANGLE RATES 

• • • 

It is convenient to calculate the Euler angle rates <() , e , and from 
the perspective of an inertial coordinate system. The mean-of-1950 (M50) 
inertial coordinate system is defined and displayed in a figure in Appendix A 
of Reference 1. Let > §2 » J .3 vectors in the directions 

of the M50 X , y , and z axes, respectively. Let 
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A = (e^- • £j) = (ij . 12 . 13 )^ 


be the orthognal transformation from M50 to body-axis coordinates. It follows 
that the ^^'s are the body axes £3 expressed in M50 

coordinates. Let 

B = (i, • Ij) ■ 


be the orthogonal transformation from M50 to topodetic coordinates. The 
£j's are the topodetic axes ^ , and expressed in M50 

coordinates. The matrices A and B are already computed within the OPIP 
program. 

The body angular velocity w expressed in body-axis coordinates is 
supplied as an input to the OPIP program. Since the body axes are fixed in 
the Orbiter, the angular velocity of the body axes relative to an inertial 
frame is the same as the angular velocity of the body. Let 

^ = A^ w 

be the body angular velocity expressed in M50 coordinates. 

Define ^ as the angular velocity of the topodetic axes as seen from the 
M50 inertial frame and as expressed in M50 coordinates. The body angular 
velocity as seen in the topodetic frame is found by removing the rotational 
motion of the topodetic frame from the inertial body angular velocity. 
Therefore 


U = W - 

is the body angular velocity as seen in the topodetic frame and as expressed 
in M50 coordinates. The next step is to compute . 

^ can be decomposed into rotational motion about f^ due to northward 
motion of the Orbiter and into rotational motions about f and f* due to 
eastward motion of the Orbiter. Let ^ and ^ be respectively the inertial 
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position and velocity of the Orbiter in M50 coordinates. An approximate value 
for ^ is 

do X R) • V di 1) • V 

^ = __^ c +_i p 

1I2 * .£1 lii III " 1^1 

(F, X R) • V 

II3.XRI |R| 

The first term above is approximate because it neglects the additional 
rotation about the increase in flattening of the Fischer ellipsoid 

with motion toward the pole. This effect is believed to be negligible. 

All rotational motion of the Orbiter as seen in the topodetic frame is 

due entirely to the Orbiter's angular velocity as seen in the topodetic 

• ^ • 

frame. Therefore, the topodetic Euler angle rates , 0 , and (j> are the 
components of body angular velocity along the instantaneous axes of rotation 
used for the topodetic Euler angles , 0 , and <|) , respectively. The next 
step is to compute the axes of rotation which will be denoted by 
hj^ , h^ , and h^ . 

The yaw is about 



The pitch is about after rotation about h^ by i|< . Since the columns 
of I (i|)) are the rotated axes in topodetic coordinates, the second column of 
gives > the rotated Jf 2 » io topodetic coordinates: 

h ^2 = (“Sini|() fj + (cos»|>) . 

The roll is about f after rotation about h, by ijj followed by rotation 

—1 —1 ^ 

about h ^2 by 0 • By a similar argument, the first column of [J( 0 ) I{i|^)] 
gives h^ , the rotated f^ , in topodetic coordinates: 


